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CHAPTER  1 
INTRODUCTION 

1.1  A  NEED  FOR  A  NEW  APPROACH  TO  ELECTROMAGNETIC 

HYPERVELOCITY  LAUNCHERS  AND  THEIR  POWER  SUPPLIES 

Electromagnetic  propulsion,  in  the  sense  of  hypervelodty  launchers 
has  been  hailed  at  the  beginning  of  this  decade  as  "an  emerging  technology 
with  the  potential  of  having  a  major  impact  in  every  aspect  of  our  lives. "[1] 

Contrary  to  the  general  expectations,  however,  the  initial  hopes  have 
not  materialized  yet,  the  projectile  velocity  records  remaining  essentially 
unbroken:  for  railguns  5.7  km /sec  obtained  by  Richard  Marshal  and 
coworkers  at  the  Australian  National  University  at  Canberra  (1970)  [2],  and  for 
coilguns,  4.9  km/sec  obtained  by  V.N.  Bondoletov  and  coworkers  in  1975  in 
the  Soviet  Union.  [3] 

A  relative  improvement  has  occurred  in  the  kinetic  energy  stored  in 
the  moving  member  of  the  electromagnetic  accelerator  [4,5],  but  at  markedly 
lower  velocities.  If  this  is  the  direction  of  improvement,  then  energies  stored 
in  the  electric  trains  built  with  linear  electric  motors  in  Japan  and  Western 
Europe  will  surpass  any  such  records. 

One  of  the  reasons  why  the  current  state  of  the  art  in  electromagnetic 
launchers  and  associated  power  supplies  technology  has  not  reached  the 
expectations  is  the  fact  that  most  of  the  activity  in  this  domain  has  been 
focused  on  the  easiest  solutions  of  simple  railgun  accelerators  (requiring  the 
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fabrication  technology  available  in  virtually  any  laboratory  or  shop)  which 
were  manufactured  predominantly  by  laboratories  which  had  power  supplies 
already  built  for  other  programs  (homopolar  generators  in  conjunction  with 
an  inductor  and  an  opening  switch,  or  high  performance  capacitors). 

The  choice  for  homopolar  generators  is  also  questionable  since  they  are 
known  for  their  low  voltage  and  low  power  density.  A  recent  study  (under 
auspices  of  D.C.  Hardison  Panel)  has  shown  that  homopolar  machines  are  not 
well  suited  for  electromagnetic  launching  applications.  (Elliott  [6]).  The  actual 
power  supply  consists  not  only  of  the  homopolar  machines,  but  includes 
heavy  and  voluminous  inductors  in  which  magnetic  energy  is  transiently 
stored  and  explosive  opening  switches  which  actually  provide  the  high 
voltage  necessary  in  order  to  drive  high  currents  into  railguns  (though  at  a 
very  low  efficiency  compared  to  alternative  capacitor  schemes).  For  Richard 
Marshall  and  coworkers  at  Canberra  [2],  the  homopolar  generator  was  existent 
and  the  only  choice. 

A  smaller  effort  was  directed  toward  heteropolar  structures  used  as 
accelerators  in  the  form  of  induction  coilguns  and  synchronous  commutated 
coilguns.  Concepts  and  preliminary  designs  have  been  performed  for  such 
alternative  structures,  but  no  notable  experiment  has  been  built  or  tested  yet. 
The  only  ones  tested  have  exhibited  less  than  mediocre  results  [7], 

Thus,  there  is  a  need  for  a  scientific  look,  from  a  fundamental  point  of 
view,  at  the  general  problem  of  electromagnetic  accelerators  and  their  power 
supplies. 
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The  analogy  with  the  electrical  machine  counterparts  -  rotating  or  even 
linear  -  has  multiple  limitations  and  can  be  counterproductive  or  even 
misleading,  due  to  the  highly  transient  and  short  duration  operation  of 
accelerators  as  compared  with  the  mostly  steady-state  performance  of  regular 
electrical  machines. 

a.  The  first  basic  feature  of  all  electromagnetic  launchers  is  the 
requirement  of  complete  integration  with  their  power  supplies. 
At  each  moment  of  operation,  the  stator  barrel  of  the  accelerator, 
the  projectile,  and  the  power  supply  must  form  a  carefully  tuned 
system,  capable  of  achieving  a  uniform  and  efficient 
electromechanical  energy  conversion  with  no  contacts  or 
feedback  control.  The  tuning  of  the  system  varies  during  the 
launch  since  the  parameters  of  the  system  must  change  in  order 
to  achieve  an  almost  constant  acceleration  and  uniformly 
distributed  body  forces  and  heating  in  the  armature  as  well  as  in 
the  barrel  as  a  means  of  surviving  the  mechanical  and  thermal 
stresses  [8]. 

b.  Another  basic  feature  of  electromagnetic  launchers  and  their 
power  supplies  is  the  variability  of  their  structure  [9].  The 
classification  patterns  and  the  characterization  following  a 
certain  classification  in  the  domain  of  classical  electrical 
machines  do  not  translate  identically  to  the  transient  nature  of 
accelerators.  For  example,  the  synchronous  and  asynchronous 
regimes  trespass  into  each  other  during  operation  and  the 
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customary  methods  of  analysis  for  one  or  another  regime  cannot 
be  used;  in  commutated  launchers,  each  commutation  occurs  at 
a  different  velocity  and  position  and  at  different  values  of  the 
currents  flowing  inside  of  the  projectile  and,  thus,  all  the 
classical  patterns  of  analytical  treatment  of  electromechanical 
devices  become  obsolete  and  not  directly  applicable. 

c.  The  phenomena  of  flux  compression  existing  only  in  some 
transient  regime  in  classical  electrical  machines  play  a 
fundamental  role  in  both  launchers  and  electromechanical 
power  supplies  of  high  performance  and  efficiency. 

d.  In  order  to  completely  define  and  characterize  the  highly 
transient  operation  of  both  accelerator  and  power  supply,  finite 
element,  three-dimensional,  transient,  electromagnetic 
computer  codes  are  necessary.  The  survival  of  the  projectile,  the 
compactness  of  the  launcher,  and  complex  problems  of  shielding 
of  the  electronics  in  the  payload  for  "smart"  projectile  require 
the  exact  knowledge  in  time  and  space,  at  any  point  in  the 
launcher  barrel,  at  any  point  in  the  projectile,  and  at  any  relative 
position  between  them  of  the  electromagnetic  fields.  From  there 
thermal  fields  and  ensuing  mechanical  stresses  can  be  calculated 
readily  transferring  data  to  similar  thermal  and  stress  finite 
element  codes  [10]. 

What  is  the  most  efficient  path  toward  hypervelocitv  electromagnetic 
launchers  and  what  are  the  power  supplies  which  best  conform  to  such 
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accelerators?  In  order  to  answer  this  question  we  will  concentrate  on  the 
"promising”  heteropolar  structures  and  will  relegate  to  other  treatments  and 
publications  about  generalized  railguns  in  which  some  high  impedance 
variants  may  bring  some  hope  of  improvement  in  the  performance  of 
homopolar  structures  used  as  launchers.  (IEEE  Transactions  on  Plasma 
Science  No  3,  June  1989,  is  devoted  almost  exclusively  to  railguns.) 

The  Bondoletov  experiment,  obtaining  4.9  km/sec  in  one  single  stage, 
shows  the  efficiency  of  repulsion  coils,  in  attaining  hypervelocities,  for 
macroparticle.  (This  reminds  us  of  the  fact  that,  in  the  realm  of  conventional 
electrical  machines,  the  repulsion  induction  motors  are  capable  of  high 
velocities,  in  very  efficient  schemes.) 

There  are  several  definite  ideas  as  points  of  reference  on  the  most 
efficient  path  toward  hypervelocity  electromagnetic  launchers. 

a.  The  first  idea  is  a  natural  thought  derived  from  the  experiment 
[3]  mentioned  above:  to  ir crease  the  number  of  stages  from  one 
to  several,  obtaining  a  "Generalized  Bondoletov  Accelerator" 
(GBA). 

b.  In  order  to  obtain  a  uniform  force  profile  in  the  projectile  and 
the  barrel,  a  current  must  be  impressed  in  the  projectile  before 
the  acceleration  had  begun  [11].  A  liquid  nitrogen  cooled 
projectile  has  the  potential  of  maintaining  the  persistent  current 
during  the  short  time  of  the  launching.  At  liquid  nitrogen 
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temperature,  the  specific  heat  has  still  an  important  value,  while 
the  resistance  of  the  copper  decreases  more  than  six  times. 

c.  Another  fundamental  characteristic  of  electromagnetic 
hypervelocity  accelerators  is  not  only  the  fact  that  the  space 
around  the  projectile  where  the  electromechanical  power 
conversion  takes  placets  strictly  localized  to  a  small  traction  of 
the  total  length  of  the  accelerator,  but,  also,  that  the  parameters 
of  such  conversion  are  continuously  changing  from  the  breech 
toward  the  muzzle  of  the  launcher. 

Then,  using  different  power  supplies  for  different  segments  of 
the  launcher  stator  (barrel)  in  the  moments  when  the  projectile 
is  contained  in  the  respective  segment  (Driga  [11])  improves 
drastically  the  efficiency  of  the  system.  Discrete  coils  in  the  GBA 
tire  pulsed  by  dedicated  electrical  machines  at  the  proper 
intervals  of  time. 

d.  It  is  shown  by  Driga,  et.  al  in  [8]  and  [12]  that  the  traveling  wave 
for  the  induction  or  synchronous  launchers  must  accelerate  at  a 
determined  pace  with  respect  to  the  projectile.  It  is  preferable  to 
obtain  such  an  accelerated  wave  by  a  continuous  increase  in  the 
frequency  of  the  power  supply,  rather  than  modifying  the  polar 
pitch  of  the  winding.  When  pulses  of  voltage  and  current  are 
used  for  discrete  coils,  for  example,  this  requirement  is  modified 
into  decreasing  the  pulse  width  and  increasing  its  amplitude  as 
the  projectile  is  accelerated. 
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The  process  can  also  be  described  a  generalized  rising  frequency  of  the 
power  supply.  Our  notion  of  frequency  variation  (learned  from  the 
customary  electrical  machines,  e.g.  the  starting  of  synchronous  motors  using  a 
varying  frequency  power  supply  with  many  cycles)  must  be  modified  for  the 
case  when  only  one  cycle,  or  even  less,  is  present.  A  generalization  in  this 
respect  means  that  the  shape  of  the  voltage  and  current  must  vary 
parabolically  -the  parameters  of  the  parabola  being  determined  by  the 
acceleration  values  and  patterns  -  at  the  given  individual  stator  coil 
(Driga[13]). 
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CHAPTER  2 

ANALYTICAL  AND  NUMERICAL  TREATMENT  OF  THE 
ELECTROMAGNETIC  LAUNCHERS  AND  THEIR  POWER  SUPPLIES 
INCORPORATING  MAGNETIC  FLUX  COMPRESSION 

2.1  A  DISCUSSION  ON  COMPENSATION  AND  FLUX  COMPRESSION 

In  the  recent  development  of  electromagnetic  launch  technology, 
electromagnetic  launchers  of  higher  impedance  (coilguns  -  as  heteropolar 
converters,  or  multitum,  augmented  railguns)  have  not  achieved  the  same 
advances  as  railgun  launchers.  There  are  several  reasons  for  this.  One  is  the 
complexity  of  such  high  impedance  devices,  where  the  payoffs  of  high 
efficiency,  high  performance,  and,  especially,  flexibility  are  counterbalanced  by 
extreme  demands  in  computations  and  design,  requiring  sophisticated 
transient  3D  electromagnetic  computer  codes  series  connected,  or  interacting, 
with  3D  stress  and  thermal  finite  element  codes.  These  codes  require  high- 
performance  supercomputers  and  the  personnel  which  few  laboratories  have 
been  able  to  dedicate  to  this  task. 

A  second,  even  more  demanding  reason,  is  the  need  for  dedicated, 
high-voltage  power  supplies  with  capabilities  for  complex  pulse  tailoring, 
advancing  considerably  the  technology  of  electromechanical  pulsed  power 
converters. 

Recent  studies  (Elliott  [6])  have  shown  that  the  inherently  low-voltage 
homopolar  machines  are  not  well  suited  for  many  electromagnetic  launching 
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applications.  The  actual  power  supply  consists  not  only  of  the  homopolar 
machine,  but  includes  heavy  and  voluminous  inductors  in  which  magnetic 
energy  is  stored  by  discharge  of  the  homopolar  machine  and  explosive 
opening  switches  which  actually  provide  the  high  voltage  necessary  in  order 
to  drive  high  currents  into  railguns  (though  at  a  very  low  efficiency, 
compared  to  alternative  capacitor  schemes). 

This  combination  of  low-voltage  rotating  machines,  matching 
inductors,  and  opening  switches  cannot  satisfy  the  waveshape  complexity 
required  by  the  higher  impedance  coilguns,  even  in  a  condition  of  extensive 
modularity. 

A  step  ahead  is  represented  by  the  Kaman  alternator  [6]  which,  being  a 
higher  voltage  machine,  leads  to  considerable  improvement  of  efficiency 
when  used  in  a  similar  scheme,  due  to  a  more  rapid  charging  of  the  inductive 
store. 


High-voltage  capacitors  achieve  higher  efficiency  and  recovery  of  an 
important  part  of  the  discharge  energy.  Their  drawback  lies  in  the  fact  that 
the  limits  of  energy  storage  densities  in  electric  fields  are  orders  of  magnitude 
lower  than  those  for  energy  densities  possible  in  magnetic  fields.  Also,  they 
require  elaborate  power  conditioning  schemes,  even  for  simple  railguns. 
Complex  power  conditioning  systems  for  the  complex  waveshapes  required 
by  heteropolar  launchers  such  as  coilguns  will  have  an  adverse  impact  to  the 
efficiency  of  the  of  the  capacitor  schemes. 
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Compensated  pulsed  alternators  have  the  advantages  of  high  voltage, 
high  efficiency,  and  capability  of  intermediate  storage  of  energy  in  magnetic 
fields.  They  are  frequently  described  as  a  one-element  power  supply  (Figure 
1),  or  single-element  power  supply  philosophy,  in  papers  by  Driga  [8  ]  and 
Driga,  et  al.  [13],  replacing  and  modifying  complex  multielement  energy 
conversion,  storage,  and  power  conditioning. 

However,  the  devices  described  above  are  not  truly  compensated 
alternators  since  both  general  procedures  for  obtaining  the  complex  pulse 
shaping  by  topological  distributions  (mainly  harmonic  synthesis)  and  by 
dynamic  interactions  (mainly  by  nonuniformly  distributed  shields  and  other 
types  of  compensation,  by  change  of  alignment  of  different  machines 
windings,  as  well  as  the  use  of  a  pulse  transformer  incorporated  in  machine), 
represents  a  means  to  decompensate,  to  undo  the  machine  compensation. 

The  descriptions  of  these  machines  includes  flux  compression  elements 
which  supposedly  provide  "active  compensation."  Actually,  flux 
compression  is  not  compensation,  as  will  be  shown  below.  However,  rotating 
flux  compressors  have  all  the  promising  characteristics  to  become  the  primary 
choice  for  power  supplies  for  high  power,  complex  heteropolar  launchers 
(complex  coilguns)  for  military,  space,  and  industrial  applications. 

The  computational  efforts  for  the  design  of  rotating  magnetic  flux 
compressors  are  tremendous  in  order  to  determine  the  time  evolution  of 
complex  interacting  electromagnetic  fields,  the  transient  mechanical  stresses, 
and  thermal  regimes.  However,  the  payoffs  are  equally  tremendous.  These 
machines  can  make  the  difference  in  assuring  the  success  of  high-power 
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Figure  1.  One-element  power  supply  philosophy 
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heteropolar  accelerators  for  hypervelocities.  On  the  other  hand,  since  the  date 
when  the  patent  [14]  for  the  compensated  pulsed  alternator,  including  the 
rotating  flux  compressor  feature,  was  issued,  advances  in  rotating  flux 
compressors  have  been  minor  (Weldon,  et  al  [37]). 

A  frequent  misunderstanding  in  approaching  pulsed  rotating  and 
linear  electromechanical  converters  is  to  equate  the  compensation  by  eddy- 
currents  produced  in  shields  in  order  to  reduce  the  internal  impedance  - 
during  the  pulsed  discharge  (by  so-called  "passive  compensation")  with  the 
flux  compression  phenomena  (incorrectly  called  "active  compensation"). 

The  first  case  is  widely  known  and  applied  frequently  in  short-circuit 
alternators  used  for  high-power  circuit  testing  of  electrical  equipment  [38].  A 
one-pulse  compensated  pulsed  alternator,  including  a  pulse-shaping  system 
based  on  harmonic  synthesis,  was  patented  by  Kapitza  in  1927  and  used  by 
Rutherford  and  Kapitza  to  produce  high  magnetic  fields  [39,  40,  41].  The 
"Pulsator"  machine  [42]  invented  by  Rebut  and  Torossian  (1966)  represents  a 
more  sophisticated  "one-pulse"  compensated  alternator  used  for  controlled 
thermonuclear  fusion  application.  Enrico  Levi  and  Harish  Pande  introduced 
a  continuous  aluminum  shield  for  compensation  and  used  an  airgap¬ 
winding  for  their  "electromechanical  pulsers"  at  the  Polytechnic  Institute  of 
Brooklyn  (1961  and  1963)  [43]. 

The  principle  of  all  these  machines  is  simple  and  is  illustrated  (Figure 
2)  below  for  a  compensated  structure  in  which  the  two-layer  winding  is  an 
airgap  winding  and  the  shield  is  continuous  and  made  of  aluminum.  The 
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eddy-currents  appearing  in  the  shield,  delay  the  penetration  of  magnetic  flux, 
and  consequently  maintain  the  internal  impedance  of  the  system  at  a  low 
level.  Figures  3,  4,  5,  and  6  show  the  unfolding  in  time  of  the  compensation 
phenomenon  at  5,  10,  15,  and  20  psec-time  intervals. 

The  second  case  is  the  flux  compression  phenomenon  and  was 
proposed  by  the  author  as  his  part  in  a  patent  for  pulsed  alternators  [14].  It  is 
achieved  in  its  simplest  form  by  an  identical  winding  connected  in  series  with 
the  armature  winding  and  placed  on  the  other  member  of  the  rotating 
machine. 

This  arrangement  does  not  provide  true  compensation  since,  during  a  360 
electrical  degree  cycle,  the  machine  takes  all  the  possible  values  for  the 
internal  impedance  from  the  maximum  to  the  minimum.  The  misnomer 
"compensation"  comes  from  the  false  analogy  with  a  dc  compensated 
machine,  where  a  similar  arrangement  provides  compensation  due  to  the  fact 
that  the  commutator  and  brush  system  maintains  the  alignment  of  magnetic 
axes  of  the  two  windings,  in  spite  of  the  rotation  of  the  armature.  In  the 
arrangement  described  below,  the  angle  between  the  two  magnetic  axes  is 
continuously  varying,  covering  all  the  possible  values  for  the  internal 
inductance  (considering  that  such  a  machine  is  "compensated"  is  analogous 
to  stating  that  a  broken  watch  shows  the  perfect  time  once  every  12  hours). 

The  essential  difference  between  "compensation"  and  "flux  compression"  is 
related  to  the  electromechanical  power  conversion  feature.  In  the  second 
case,  the  system  of  two  windings  linking  a  magnetic  flux  in  the  conditions  of 
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Figure  3.  Magnetic  flux  at  t  =  5  jisec 
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Figure  3a.  Vector  potential  values  at  t  =  5  nsec 


Contour  number  r  X  A  First  point  found  in  El moment  numbered 

r  z 


1 

1 . 71027E-08 

5.  5322 IE— 03 

-4.  00000E— 03 

77 

2 

3.  42054E-08 

6.  60000E— 03 

-3. 64420E— 03 

79 

3 

5.  13080E-08 

6.  60000E— 03 

-2. 38445E— 03 

79 

4 

6.  84107E-08 

7.  20000E— 03 

-2. 35157E— 03 

80 

5 

8.  55134E-0S 

7.  2000QE— 03 

-1. 78960E— 03 

109 

6 

1. 02616E-07 

7.  20000E— 03 

-1. 46341E-03 

138 

7 

1.  19719E— 07 

6.  60000E— 03 

-1 . 07205E-03 

166 

8 

1. 36821E-07 

7.  20000E-03 

-9. 58662E-04 

167 

9 

1. 53924E-07 

4.  13322E— 02 

-5. P1437E-04 

188 

10 

1. 71027E-07 

4.  46667E— 03 

-5. 63228E-04 

189 

1  1 

1. 8B129E -07 

3.  33333E— 03 

-5.  68966E— 04 

192 

12 

2. 05232E— 07 

4.  13333E-03 

-4.  2B142E-04 

217 

13 

2. 22335E-07 

4.  13333E-03 

-4. 62129E— 04 

217 

14 

2. 39438E-07 

4.  46667E— 03 

-4. 59439E-04 

218 

15 

2. 56540E— 07 

5. 33333E— 03 

-4. 71676E-04 

221 

16 

2. 72642E—07 

5. 66667E— 03 

-4.  5 1 507E— 04 

222 

17 

2. 90746E-07 

4. 46667E— 03 

-3. 78475E-04 

247 

18 

3. 07848E— 07 

4.  46667E— 03 

—3.  5 1094E— 04 

247 

19 

3. 2493  IE-07 

3.  33333E— 03 

-3.  74704E— 04 

230 

20 

3. 42054E-07 

4. 46667E— 02 

-2.  76764E-04 

276 

21 

3. 59156E— 07 

5. 00000E— 03 

—2. 52543E— 04 

278 

22 

3. 76259E— 07 

5. 33333E— 03 

—2. 79454E-04 

279 

23 

3. 93362E-07 

5. 00000E— 03 

-1. 53050E-04 

307 

24 

4. 10464E— 07 

5. 33333E— 03 

-1. S2612E— 04 

308 

25 

4. 27567E— 07 

5. 00000E— 02 

-2.  19852E—05 

336 

26 

4. 44670E-07 

5. 32333E— 03 

-3. 8 196 IE-05 

337 

27 

4. 61772E-07 

5. 00000E— 03 

1 . 20807E-04 

365 

28 

4.  78875E-07 

5. 33323E— 03 

9. 75367E-05 

366 

29 

4. 95978E-07 

5. 32322E— 03 

1. 66569E-04 

395 

30 

5.  13080E-07 

5. 33333E— 03 

3. 50000E-04 

424 
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Figure  4a.  Magnetic  vector  potential  values  at  t  =  10  psec 


Contour  number  r  X  A  First  point  found  in  Eimement  numbered 

r  2 


1 

3. 61410E— 08 

5.  25737E— 03  00000E-03 

76 

2 

7.  22820E-O8 

6.  19910E— 03  —3. 35292E— 03 

79 

3 

1 . 08423E-07 

6.  60000E— 03  —2. 8 406 4E— 03 

79 

4 

1 . 44564E— 07 

7.  20000E— 03  -2.  6658BE-03 

80 

5 

1 . 80705E— 07 

7.  80000E-03  —2. 22314E— 03 

81 

6 

2. 16846E-07 

7,  80000E— 03  — 1 . 71033E-03 

110 

7 

2. 52987E-07 

7.  20000E— 03  — 1 . 32633E—03 

138 

8 

2. 8912BE— 07 

6.  60000E— 03  -1.03483E-03 

166 

9 

3. 25269E— 07 

3.  60000E— 03  —5.  69286E— 04 

186 

10 

3. 61410E— 07 

4.  13323E— 03  -6. 04E73E-04 

188 

11 

3.  9755 IE-07 

4.  13333E-03  -5.56776E-04 

188 

12 

4. 33692E-07 

4.  80000E— 03  -5.53115E-04 

190 

13 

4. 69833E-07 

5.  33333E— 03  -5.59B44E-04 

192 

14 

5. 05974E-07 

4.  13333E—03  -4. 69568E-04 

217 

15 

5. 421 15E—07 

4.  46667E— 03  68438E-04 

218 

16 

5.  78256E— 07 

5.00000E-03  —4. 52977E— 04 

220 

17 

6. 14397E-07 

5.  33333E— 03  -4.  60063E-04 

221 

18 

6. 50538E— 07 

4.  46667E— 03  -3.81746E-04 

247 

19 

6. 86679E-07 

4.  46667E— 03  -3. 52069E-04 

247 

20 

7. 22820E-07 

5.  33333E— 03  -3.79866E-04 

250 

21 

7. 5896 IE-07 

5.  33333E—03  -3. 52675E-04 

250 

22 

7.  95102E-O7 

5.  33333E— 03  -3. 15237E-04 

279 

23 

8.  31242E-07 

5.33333E-03  -2.76681E-04 

279 

24 

8.  67383E-07 

5.  33333E— 03  -2.29912E-04 

308 

25 

9. 03524E-07 

5.33333E-03  -1. 64691E-04 

308 

26 

9. 39665E-07 

5.  33233E— 03  -8.  14006E-05 

337 

27 

9.  75806E-07 

4. 80000E-03  1 . 40253E— 04 

364 

28 

1. 01195E-O6 

5. 33333E— 03  7. 80595E-05 

366 

29 

1 . 04809E-06 

5.  33333E— 03  1.50401E-04 

395 

30 

1 . 08423E-06 

5.  33333E— 03  2. 50000E-04 

395 

20 


Figure  5a.  Magnetic  vector  potential  values  at  t  =  15  gsec 


Contour  number  r  X  A 


First  point  found  in  Elmement  numbered 
r  2 


1 

5. 65B45E— 08 

5.  00353E— 03 

-4. 00000E— 03 

76 

2 

1 .  13169E-07 

6. 00090E— 03 

-3. 0 1 95  IE— 03 

78 

3 

1 . 69753E— 07 

£.  60000E— 03 

-3. 1 4876E— 03 

79 

4 

2. 26338E-07 

6. 60000E-03 

-2. 1 9407E— 03 

79 

5 

2.  82922E-07 

7.  20000E— 03 

-2. 16063E-C3 

80 

6 

3. 39507E-07 

7.  20000E—03 

-1. 74285E-03 

109 

7 

3. 96091E— 07 

7.  2C000E— 03 

-1. 44E38E— 03 

138 

8 

4. 52676E-07 

6. 60000E— 03 

-1. 1 2952E— 03 

166 

9 

5. 09260E— 07 

6. 60000E— 03 

-9. 75743E— 04 

166 

10 

5. 65845E— 07 

3.  60000E— 03 

-5. 65B48E— 04 

186 

1  1 

6. 22429E-07 

4.  13333E— 03 

-6. 0751  IE-04 

188 

12 

6.  79014E— 07 

4.  13333E— 03 

-5. 60372E-04 

188 

13 

7. 35598E-07 

4. 46667E— 03 

—5. 52470E— 04 

189 

14 

7. 921B3E-07 

5.  33333E— 03 

-5. 70756E— 04 

192 

15 

8. 48767E—07 

4.  13333E— 03 

-4. 6679  IE-04 

217 

16 

9. 0535 IE— 07 

4.  46667E-03 

-4. 6751  IE-04 

218 

17 

9. 61936E— 07 

5. 00000E— 03 

-4. 53777E-04 

220 

18 

1. 01852E-06 

5. 33333E— 03 

-4. 60876E— 04 

221 

19 

1. 07510E—06 

4. 46667E— 03 

-3. 74902E— 04 

247 

20 

1.  13169E-06 

5.  00000E— 03 

-3. 57630E-04 

249 

21 

1.  18827E— 06 

5 . 33333E— 03 

-3. 75 1 48E— 04 

250 

22 

1. 24486E-06 

4.  46667E-03 

-2. 53  1 43E-04 

276 

23 

1. 30 1 44E-06 

5.  33333E— 03 

-3. 02839E-04 

279 

24 

1. 35803E— 06 

5. 33333E— 03 

—2 . 6 1 292E— 04 

279 

25 

1. 4 1461E— 06 

5.  33333E— 03 

-1. 97859E-04 

308 

26 

1. 47120E-06 

5.  00000E  -03 

-9. 46264E-06 

336 

27 

1 . 52778E-06 

P.  333 j j 

89386E-05 

337 

28 

1. 58437E-06 

5. 33333E— 03 

6. 37434E-05 

366 

29 

1. 64095E-06 

5.  33333E— 03 

1. 42674E— 04 

366 

30 

1. 69753E-06 

5.  33333E— 03 

2. 50000E— 04 

395 

21 


Figure  6.  Magnetic  flux  at  t  -=  20  (isec 
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Figure  6a.  Magnetic  vector  potential  value  at  t  =  20  psec 


Contour  nunbor  r  X  A  first  point  found  in  Eiowawnt  numbered 

r  x 


1 

8.  72344E— 08 

4.  49G38E— 03 

-4. OOOOOE-03 

74 

2 

1 . 74469E-07 

5.  666C7E— 03 

-2. 61419E-03 

77 

3 

2.  617Q3E— 07 

6. 00000E-03 

-2. 27288E-03 

78 

4 

3. 48937E— 07 

6 . 60000E— 03 

-2. B3417E-03 

79 

5 

4.  36172E— 07 

6. 60000E— 03 

-2. C5208E— 03 

79 

6 

5.  23406E— 07 

7.  80000E-03 

-2. 20496E-O3 

81 

7 

6. 10640E— 07 

7.  20000E— 03 

—1 . 7C5S3E— 03 

109 

8 

6.  97875E— 07 

C. 60000E-G3 

-1. 3418SE-03 

137 

9 

7. 85109E— 07 

Z. 60000E-03 

-7. 62785E-04 

157 

10 

8.  72344E-07 

4.  12323E— 03 

-7. 90E58E— 04 

159 

11 

9.59578E-07 

4.  4C6C7E— 03 

—7. 71904E—04 

160 

12 

1. 04681E-06 

5. 0Q00QE-03 

—7. 55330E— 04 

1C2 

13 

1. 13405E— 06 

3.  60000E— 03 

-5. 72009E-04 

186 

14 

1 . 22128E—06 

3.  80000E— 03 

—5.  6 1 852E-04 

187 

15 

1 . 30852E-06 

4.  13333E— 03 

-5. 82448E-04 

188 

16 

1 . 39575E-06 

4.  46667E— 03 

-5. 82540E— 04 

189 

17 

1 . 48298E— 06 

4.  80000E-03 

-5.  62034E-O4 

190 

18 

1 . 57022E-O6 

5. 23333E-03 

-5. 6C542E-04 

192 

19 

1 . 65745E-06 

4. 46667E-03 

-4. 76029E-O4 

218 

20 

1 . 74469E-06 

4. 80000E— 03 

-4. 51 181E—04 

219 

21 

1. 83192E-06 

5. 33333E-03 

-4. 68557E— 04 

221 

22 

1. 91916E-06 

4. 46667E— 03 

-3.  59542E|-04 

247 

23 

2. 00639E-06 

5.  33333E— 03 

-3. 96905E— 04 

250 

24 

2. 09362E-06 

5.  33333E— 03 

-3. 5973  IE-04 

250 

25 

2. 18086E— 06 

5.  33333E— 03 

-3. Q9463E— 04 

279 

26 

2. 26809E— 06 

5.  33333E— 03 

-2. 54552E— 04 

279 

27 

2. 35533E— 06 

5.  33333E— 03 

-1 . 5 5 26 7E— 04 

308 

28 

2. 44256E-06 

5.  33333E— 03 

-6.  47236E— 06 

337 

29 

2. 52980E-06 

5.  33333E—03 

1. 06247E-O4 

366 

30 

2. 6 1703E-O6 

5.  33333E— 03 

2. 50000E-04 

395 

23 


a  large  inductance  and  consequently  small  current  is  compressed  by  external 
forces,  converting  mechanical  (kinetic)  power  into  electrical  power  -  in 
condition  of  limited  flux  losses  -  until  it  reaches  the  position  of  minimum 
inductance  in  the  presence  of  large  currents. 


2.2  SIMPLIFIED  FLUX  COMPRESSION  CIRCUIT  ANALYSIS  AND 

ELECTROMECHANICAL  POWER  CONVERSION  CONSIDERATIONS 


The  classical  treatment  of  the  simplest  flux  compression  system 
involves  a  simple  R,  L  circuit,  in  which  the  flux  to  be  compressed  is  /stored 
initially  in  the  inductance  by  means  of  an  external  power  supply  which  in 
Figure  7  is  represented  as  a  current  source  which  is  already  short-circuited  at 
f=0  when  the  compression  begins.  In  principle,  the  flux  is  compressed  by 
rapidly  modifying  the  relative  position  of  an  arrangement  of  conductors  [15], 
forcing  the  magnetic  flux  to  conform  to  a  continuously  decreasing  inductance 
Lit)  .  Since  the  magnetic  energy  (as  well  as  the  magnetic  coenergy)  in  the 
system,  Wm,  can  be  expressed  as  a  function  of  the  instantaneous  current,  1(A), 
and  the  magnetic  flux,  as  state  variables: 


i  1  vt/2 

Wm=-L(t)I2 
n  2  2  L{t) 


(2.1) 


it  can  be  seen  that  in  condition  of  constant  flux,  the  magnetic  energy  stored  in 
the  system  and  partly  delivered  to  the  load  is  continuously  increasing  until 
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the  minimum  inductance  position  is  reached.  It  is  obvious  that  the  increase 
in  electromagnetic  energy  occurs  by  conversion  of  the  mechanical  energy 
stored  as  kinetic  energy  in  the  moving  conductors.  However,  the  flux  does 
not  remain  constant,  being  affected  by  losses  into  the  circuit  resistance,  R. 

The  differential  equation  for  the  circuit  in  terms  of  magnetic  flux  is: 


d(U) 

dt 


+  y(U)  =  0, 

L 


(2.2) 


where  the  magnetic  flux  V  has  been  expressed  in  terms  of  inductance: 


V  =  U(Wb) 


(2.3) 


The  solution  of  the  differential  equation  (2)  is: 


U  =  L0I0  exp 


(2.4) 


where 


L 


total  inductance  of  the  circuit  at  time  t,  (H) 
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1  =  current  through  the  circuit  at  time  t,  (A) 

Lq  =  initial  total  inductance  at  time  t  =  0,  (H) 

Iq  =  initial  current  (time  t=0),  U) 

Rl  =  the  load  resistance  which  somewhat 

increased  to  take  into  account  the  resistance 
of  the  rest  of  the  circuit  [fl  J 


As  the  total  inductance  is  decreasing,  the  current  is  increasing: 


(2.5) 


In  order  to  obtain  the  instantaneous  power  balance,  we  multiply  (2.2)  by  I,  and 
by  rewriting  the  first  term,  we  get: 


2  dx  dt 


KM 


(2.6) 


Considering  that  the  flux  compression  is  obtained  by  a  translational  motion 
in  the  direction  of  the  generalized  coordinate,  x„  we  obtain: 


1  j 2  dL  dx' 

2  dx__dt_ 


=  RI2  +  — 
dt 


(2.7) 


where  the  two  pairs  of  brackets  in  the  left  hand  member  contain  the 
generalized  force  on  the  x  direction  fand  the  velocity  in  the  same  direction, 
u  Respectively: 
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-F*u=Rl2 


dW 

+ - 

dt 


(2.7’) 


The  mechanical  power  input  necessary  to  decrease  the  inductance  of 
the  circuit  is  spent  in  increasing  the  magnetic  energy  stored  in  the  system  and 
is  partially  lost  in  Joule  heat. 


If,  as  in  rotating  machinery,  a  rotational  motion  is  responsible  for  the 
flux  compression  by  decreasing  the  inductance  of  the  circuit,  relation  (2.7) 
becomes: 


'da' 

.2 1  da. 

• 

.  dt . 

=  RI2+  — 
dt 


(■ 


- U 2 
2 


(2.8) 


where,  in  the  left  side  we  have  the  torque  T(N*m)  as  generalized  force  in  the 
a  direction  and  the  angular  velocity  5J  rad/s)  as  the  derivative  of  the  angular 
displacement  arad  with  respect  to  time.  Then  the  instantaneous  power 
balance  becomes: 


-T*(d  =  RI 2 


+ 


dWn 

dt 


(2.9) 


where  Wm  is  the  magnetic  energy  stored  in  the  circuit  at  that  particular 
moment. 


However,  pulsed  electrical  generators  incorporating  flux  compression 
features  have  the  advantage  of  producing  new  flux  during  compression 
which  is  continuously  integrated  and  compressed,  in  addition  to  the  initial 
flux.  Then  the  equation  (2.2)  becomes 
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|[ZJ]+|[L/]  =  /«> 


(2.2') 


The  simplicity  of  the  equation  is  deceiving,  f(t)  being  a  complicated 
function,  since  the  excitation  flux  and  the  angular  velocity,  (o/  are  variable 
and  nonlinear  during  the  compression  and  the  differential  equation  must  be 
linked  to  the  instantaneous  balance  of  power,  and  to  the  diffusion  equation 
for  the  electromagnetic  fields  in  3D  structures  requiring  elaborate  numerical 
codes,  in  condition  of  relative  motion. 


In  equations  (2.7-2.9)  the  converted  power  appearing  in  the  left  side 
member  has  a  negative  sign.  In  electromagnetic  converters,  this  sign  shows 
the  direction  of  energy  flow.  This  case  indicates  the  generator  -  flux 
compressor  regime  in  the  linear  topology  variant,  equations  (2.7)  and  (2 .7’) 
and  in  the  rotating  variant,  equations  (2.8)  and  (2.9). 

The  mechanical  energy  comes  from  the  kinetic  energy  usually  stored  in 
the  mass  of  moving  conductors  achieving  compression.  In  turn,  the  kinetic 
energy  is  produced  by  some  form  of  primary  energy  coming  from  an 
explosive  for  the  common  flux  compressors  designed  to  achieve  ultra-high 
fields  or  from  a  turbine,  in  the  case  of  electromechanical  power  supplies  for 
launchers  (e.g.  compulsator)  [14]. 

Conversely,  for  electromagnetic  accelerators,  the  terms  describing  the 
converted  power  and  the  magnetic  energy  stored  in  the  circuit  change  signs  in 
order  to  show  the  change  in  direction  of  energy  flow,  from  magnetic  energy  to 
an  increase  in  the  kinetic  energy  of  the  projectile: 


29 


F*u=RI 2 


dWm 

dt 


(2.7") 


and 


T  •  co  =R12 


dWm 

dt 


(2.9') 


The  equations  (2.1)  to  (2.9)  are  global  or  integral  relations  describing  the 
whole  magnetic  flux  compression  (generators)  or  accelerator  circuits.  Point 
(or  local)  relations  may  be  used  for  a  local  and  distributed  characterization  of 
the  power  conversion. 

Following  the  classical  treatment  of  Levi  and  Panzer  [16]  Nasar  [17],  or 
White  and  Woodson  [18],  the  unit  volume  element  of  the  conductor  is 
considered  the  building  block  of  the  converter.  By  establishing  the  different 
specific  power  densities  and  by  multiplying  them  by  the  elementary  volume 
dw,  the  electromechanical  energy  conversion  process  is  characterized  and 
defined  locally,  in  each  point  of  the  converter. 


The  power  balance  relations  for  the  unit  volume  of  conductor  [16]  are: 
C  «•  — =  (JxJ?)»S  +  /  •« 

S  dt  (2.10) 

J2  —  -  —  - 

—  =  E*J  +(uxB)*J 

O  (2.11) 


We  denote  the  power  per  unit  volume  or  specific  power  with  lower 
case  p(w/m 3).  Then,  =  specific  kinetic  power,  Pco*>=  specific  converted 
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power  transformed  into  heat  by  Joule  effect.  Then: 


r-du  d  (1  y 

ptim=fu  —  =  —  —Cu 
fa"  dt  dt  1 2 


Pco„=OxB)u 

pn<ck=L*“ 
p„  =  E*7 
J2 

Pbss- - 


(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 


leading  to  the  equivalent  circuits  per  unit  volume  of  conductor  of  a  magnetic 
flux  compressor  (generator)  (Figure  8)  and  of  an  electromagnetic  accelerator 
(Figure  9). 

What  must  be  the  law  of  variation  for  velocity  u  ?  One  of  the 
optimum  situations  in  electromagnetic  accelerator  design  is  reached  when 
the  forces  are  uniformly  distributed  over  the  projectile  -  assuring  a  uniform 
distribution  of  stresses  -  and  leading  to  a  constant  acceleration  during  the 
launch  time.  In  [8]  (Driga,  et.al.),  the  means  to  achieve  an  accelerated 
travelling  wave  following  almost  exactly  the  motion  law  of  the  projectile  by 
using  a  continuously  rising  frequency  or  a  variable  pitch  winding  or  a 
combination  of  both  was  described. 


As  shown  in  Driga  [9],  a  rising  frequency  generator,  while  possible  in 
different  variants  in  a  pulsed  mode,  is  a  very  demanding  and  complex 
electrical  machine  which  works  properly  only  for  a  limited  span  of  its 
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circumference.  But  even  then,  such  a  power  supply  would  not  create  the 
purely  sinusoidal  travelling  wave  acting  uniformly  on  the  currents  in  the 
projectile  -  ideal  condition  which  is  unattainable  in  practice.  The  reasons  for 
that  [19]  are  related  to  the  harmonics  arising  due  to  the  finite  length  of  the 
barrel  wave  packet,  its  lack  of  perfect  smoothness,  and  the  finite  length  of  the 
projectile  whose  leading  and  lagging  ends  must  slip  ahead  and  behind  the 
wave  packet,  only  the  middle  of  it  travelling  at  the  speed  of  the  accelerated 
wave. 
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Figure  8  Electromechanical  power  conversion  in  a  unit  volume  of  an 
electromagnetic  flux  compressor 


/ 


Figure  9.  Electromechanical  power  conversion  in  a  unit  volume  of  an 
electromagnetic  launcher 
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2.3  ANALYTICAL  AND  NUMERICAL  TOOLS  FOR  ANALYSIS  AND 

SYNTHESIS  OF  ELECTROMAGNETIC  LAUNCHERS  AND  POWER 

SUPPILES 

Then,  what  are  the  right  tools  for  the  analysis  and  synthesis  of  power 
supplies  and  electromagnetic  accelerators?  Obviously,  the  rapid  evolution  of 
the  magnetic  flux  distributions  and  currents  accompanying  them  needs  a 
global  characterization;  while  locally,  the  evaluation  of  the  flux  densities  and 
current  densities  describe  the  electromagnetic  state  at  each  point. 

It  is  widely  assumed  that  the  next  breakthroughs  and  innovations  in 
hypervelocity  accelerator  technologies  will  come  from  advances  in  analysis 
and  modelling,  controlling  electromagnetic  diffusion  and  fast  transient 
phenomena  in  the  projectile  and  launcher  barrel,  and  in  the  power  supplies. 
This  leads  us  to  consider  the  finite  element  method  as  the  main  tool  for 
analyzing  the  transient  electromagnetic  phenomena  during  the  magnetic  flux 
compression  -  expansion. 

Later  in  Chapter  3  the  above  problem,  involving  a  3-D  f.e.m. 
formulation  using  a  Galerkin  technique,  is  described  in  detail. 

Driga,  et.  al.,  [10]  and  Pillsbury  [24]  describe  a  2-D  axisymmetric  f.e.m. 
formulation  of  a  transient  electromagnetic  diffusion  in  fast  discharging 
homopolar  machines,  also  used  to  illustrate  simple  magnetic  compensation 
arrangements. 
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Magnetic  flux  compression  -  expansion  devices  require  the  knowledge 
of  3-D  flux  distributions,  and  since  the  magnetic  flux  and,  consequently, 
magnetic  stress  concentrations  occur  at  connections  where  the  field  is  three- 
dimensional,  simplifications  to  two-dimensional  structures  are  no  longer 
acceptable. 

The  computational  efforts  for  a  three-dimensional  f.e.m. 
electromagnetic  computer  code  are  very  large  since  such  a  transient  code 
solves  the  problems  for  all  three  components  of  the  vector  potential 
(Ax,Ay,At). 

By  comparison,  the  similar  code  for  two-dimensional  axisymmetric 
structures  mentioned  above  solves  the  problem  for  only  one  component  of 
the  magnetic  vector  potential  (A® ). 

One  of  the  most  important  features,  important  not  only  for  numerical 
reasons,  but  having  also  deep  theoretical  implications,  is  the  selection  of  the 
Coulomb's  gauge  for  the  magnetic  potential: 

divA  =  0  (2.12) 

This  choice,  beside  completely  defining  the  field  of  vectors  A ,  makes  the 

statement  that  the  quasistationary  assumption  dt  holds.  Also,  that 
locally: 


divJ  =  0 


(2.13) 
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is  valid  for  all  the  magnetic  flux  compression  -  expansion  devices  (or, 
globally,  in  circuit  terms,  the  first  Kirchhoffs  law  applies).  Equation  (2.12)  is 
immediately  derived  from  (2.13),  taking  into  account  that  the  product  R2\A\ 
remains  finite  at  large  distances  ( R  — > «).  Coulomb's  gauge  choice,  required 
point- wise,  is  justified  by  Driga,  et.  al.  [10]  and  Pillsbury  [24]  as  simplifying  the 
governing  equations  for  the  2-D  formulation,  and  its  application.  In  that 
particular  code,  was  imposed  by  the  use  of  Lagrange  multipliers.  However,  in 
the  axisymmetric  geometry,  the  condition  follows  naturally,  since  prescribing 
a  current,  the  entire  circuit  is  implied,  and  condition  (2.13)  is  respected.  Also, 
in  such  conditions,  the  Lagrange  multipliers  may  be  replaced  equivalently  by 
some  appropriate  penalty  function. 

In  the  3-D  finite  element  code,  the  situation  changes  completely.  In 
such  geometries,  the  usual  forcing  function  is  prescribed  brick  by  brick,  and 
small  discrepancies  introduced  by  round-off  errors,  for  example,  lead  to  the 
fact  that  Coulomb's  gauge  becomes  invalid.  Also,  the  nonidentical  three- 
dimensional  elements  joining  each  other  add  to  discrepancies  in  the 
condition  (2.12),  point-wise. 

The  tool  imposing  the  Coulomb's  gauge  constraint  globally  per 
element  and,  also,  point-wise  in  the  sense  of  distribution  along  two  adjacent 
integration  points,  is  the  use  of  Lagrange  multipliers. 

As  shown  in  Chapter  4,  the  Lagrange  multipliers  in  this  three- 
dimensional  code  represent  sheets  of  currents  flowing  through  the  faces  of 
the  three-dimensional  elements  and,  by  compensating  the  discrepancies. 
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impose  the  Coulomb's  gauge,  point- wise,  in  the  entire  volume.  The 
numerical  and  theoretical  implications  of  this  fact  represent  some  of  the 
original  contributions  of  this  work. 

The  nonuniformities  in  prescribing  the  forcing  function  and  other 
nonuniformities  introduced  by  the  nature  of  the  finite  element  method  itself 
(finite  elements  for  example...)  are  compensated  by  the  magnetic  flux 
produced  by  surface  currents  flowing  through  the  faces  of  the  volume 
elements.  Then,  in  Chapter  5,  this  compensation  method  is  applied  to  the  3- 
D  finite  element  formulation  for  solving  the  complete  system  of  Maxwell's 
equations,  under  the  constraint  divA  =  0.  The  three-dimensional,  transient 
code  is  applied  to  a  generalized  Bondoletov-type  accelerator,  having  both 
massive  conductive  elements  characterized  by  free  diffusion  and  stranded 
and  transposed  conductive  elements,  in  which  the  flow  of  current  is 
controlled.  The  magnetic  field  and  current  density  space  distributions  are 
plotted  at  successive  intervals  for  a  rectangular  stator  coil  repelling  a 
rectangular  copper  plate. 

Comparisons  are  made  with  a  similar  flux  compression  system  with 
round  coil  and  round  plate  which,  due  to  its  axisymmetric  arrangement,  can 
be  solved  with  a  two-dimensional  code.  The  difference  between  the  two 
solutions  is  theoretically  fundamental,  being  based  on  the  manner  in  which 
the  constraint  divA  =  0  is  imposed. 

In  the  latter  case,  it  follows  almost  automatically  due  to  the  symmetry 
of  the  current  (and  vector  potential)  distributions;  in  the  first  case,  it 
practically  illustrates  the  considerations  made  in  Chapter  4,  showing  how 
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efficient  are  the  limit  cases  of  sheets  of  currents  flowing  along  the  faces  of  the 
volume-finite  elements  in  establishing  the  right  (quasistationary)  distribution 
of  magnetic  flux. 

The  computational  effort  for  analyzing  and  synthesizing  magnetic  flux 
compression  -  expansion  convertors  using  using  three-dimensional  f.e.m. 
transient  codes  is  tremendous.  The  payoffs  are  equally  tremendous,  since 
only  perfectly  designed  defices  can  assure  the  success  of  hypervelocity 
accelerators,  overcoming  the  extreme  stress  levels  with  an  exact  knowledge  of 
their  transient  evolution. 

The  examples  presented  are  three-dimensional  plots  of  magnitude  and 
directions  for  magnetic  field  and  currents  density  of  a  flux  compression 
system  with  transient  field  diffusion.  The  consecutive  stages  of  the  magnetic 
flux  compression  are  shown  at  consecutive  intervals  of  5  psec,  namely  5,  10, 
15,  and  20  psec.  (During  this  time  the  inertia  causes  the  projectile  to  move 
very  little.) 

There  is  another  important  reason  why  the  two-dimensional 
axisymmetric  codes  are  ineffective  and  must  be  replaced  by  three-dimensional 
ones  (even  if  apparently  small  changes  in  structure  can  be  made  to  the  latter 
to  look  axisymmetric),  namely  the  armature  reaction.  The  armature  reaction, 
changing  in  time  during  launching  due  to  its  spatial  shift  of  the  magnetic 
fields,  adds  irrevocably  the  third  dimension  component. 

The  output  of  the  3-D  code,  in  terms  of  magnetic  vector  potential, 
permits  the  calculation  of  magnetic  fields  and  current  density  distributions. 
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electromagnetic  body  forces,  Joule  losses,  and  the  density  of  magnetic  energy 
stored  in  the  field.  Then,  all  terms  in  the  formulas  (2.12)  to  (2.16)  can  be 
calculated  and  the  electromechanical  power  conversion  characterized  in  each 
point  of  the  electromagnetic  device. 

Globally,  by  integration,  the  magnetic  flux  being  compressed  results 
almost  directly,  since  the  vector  potential,  A,  and  the  flux  density,  B ,  are 
related  by  the  Stokes'  theorem: 

<fr  =  j^B»ds  =  §A*d£ 

S'  C 

where  the  surface  ■$.  lies  on  the  curve  c. 

Similarly,  the  currents,  the  total  magnetic  energy  stored  in  the  field, 
total  Joule  losses,  and  global  forces  car.  be  calculated  from  local  field  values 
and  integrated. 

Then,  the  formulas  (2.5)  to  (2.9)  may  be  applied  characterizing  the 
power  conversion  phenomena  globally  in  the  magnetic  flux  compression  - 
expansion  devices. 

In  fact,  the  forces  of  electromagnetic  origin  (and  the  velocities  which 
result  from  the  continuous  acceleration  in  the  barrel)  are  of  main  concern  in 
obtaining  new  limits  of  performance  in  electromagnetic  hypervelocity 
accelerators.  Such  forces  can  be  calculated  in  three  different  ways. 

1.  As  the  global,  resultant  force  (a  generalized  Lagrangian  force, 
acting  along  a  generalized  coordinate); 
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2.  As  the  distribution  of  body  forces  in  the  projectile  volume  and 
in  the  rest  of  the  accelerator; 

3.  And,  finally,  as  stresses,  using  the  artificial  method  of 
electromagnetic  stress  tensor,  reducing  the  problem  of  electric 
and  magnetic  fields  of  force  to  that  of  an  elastic  continuum. 

The  three  methods  are  formally  equivalent  and  their  expressions  result 
directly  and  elegantly  due  to  the  almost  symmetrical  structure  of  Maxwell's 
equations  for  macroscopic  nonrelativistic  electrodynamics. 

The  theoretical  tools  which  transform  the  expressions  for  forces  from 
one  form  to  another  are  related  to: 

a)  The  constancy  of  magnetic  flux  as  the  constraint  used  to  assure 
the  proper  transformation  in  order  to  find  the  generalized 
Lagrangian  forces  along  a  generalized  coordinate  (global  or 
integral  characterization);  and 

b)  The  complex  and  subtle  notion  of  flux  derivative  (2.12)  as 
applied  to  electromagnetic  systems  in  motion.  Again,  using  the 
constancy  of  the  magnetic  flux  as  a  constraint,  we  obtain  all  the 
equivalent  body  forces  of  electromagnetic  origin  (local  or  point 
characterization). 

Such  results  obtained  in  Chapter  5  permit  a  direct  transition  from  the 
global  to  local  (and  reverse  from  local  to  global)  characterization  of  forces  in 
electromagnetic  accelerators  and  their  power  supplies. 


40 


CHAPTER  3 

THE  FORMULATION  OF  THE  FINITE  ELEMENT.  THREE  DIMENSIONAL 
TRANSIENT  ELECTROMAGNETIC  PROBLEM  FOR  MAGNETIC  FLUX 
COMPRESSION  DEVICES  AND  ELECTROMAGNETIC  LAUNCHERS 

3.1  INTRODUCTION 

The  finite  element  method  is  ideal  for  analyzing  the  transient 
magnetic  field  diffusion  in  the  flux  compression  devices  in  three- 
dimensional  domains.  It  can  easily  handle  discontinuous  geometrical  shapes 
and  material  discontinuities  in  solving  the  penetration  of  magnetic  flux  into 
complex  structures  of  the  flux  compression  devices. 

A  formulation  of  three-dimensional  transient  electromagnetic 
problem  for  flux  compression  is  presented  in  this  chapter.  It  solves  the 
magnetic  field  diffusion  problems  in  terms  of  all  three  components  of  the 
magnetic  vector  potential  by  comparison  with  the  similar 

numerical  codes  for  two-dimensional  axisymmetric  electromagnetic 
structures  which  solve  the  flux  penetration  problems  for  one  component  of 
the  magnetic  vector  potential 

Magnetic  flux  compression  devices  comprise,  among  various 
conductive  and  insulating  materials,  also  stranded  and  transposed 
conductors.  These  components  are  also  taken  into  account  by  the  present 
formulation  by  the  artifice  of  considering  their  conductivity  zero  and 
prescribing  their  current  densities  as  a  function  of  time,  usually  obtained 
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with  good  accuracy  using  circuit  analysis  methods,  or,  alternatively,  from  an 
iteration  procedure. 

One  of  the  important  characteristics  of  the  present  formulation  is  the 
introduction  of  superficial  currents  -  flowing  through  the  faces  of  the  three- 
dimensional  finite  elements  -  as  Lagrange  multipliers. 

The  approximation  of  the  magnetic  flux  which  is  a  perfect  continuum 
due  to  the  solenoidal  characteristics  of  the  flux  density  B,  by  the  f.e.m,  a 
discrete  method  of  solving  parabolic,  partial  differential  equations  [30],  is 
performed  by  imposing  the  Coulomb  gauge  [31],  pointwise.  This  is  based  on 
the  fact  (shown  in  Chapter  4)  that  the  magnetic  field  produced  by  a 
distribution  of  volume  currents  inside  an  element  (brick,  prism,  etc.)  can  be 
replaced  by  an  identical  magnetic  field  structure  produced  by  a  distribution  of 
surface  currents  [26],  [31]  flowing  on  different  faces  of  the  respective  elements. 
The  prescription  of  current  densities  element  by  element,  round-off  errors  in 
calculating  the  current  densities,  the  slight  mismatch  of  adjacent  or  boundary 
elements,  etc.,  will  require  a  strong  imposition  of  Coulomb’s  gauge  in  order 
to  enforce  the  solenoidal  character  of  the  current  density  everywhere. 

3.2  DOMAIN  OF  DEFINITION  AND  GENERAL  ASSUMPTIONS 

The  domain  under  consideration,  extending  eventually  to  infinity  in 
all  three  dimensions,  has  (Figure  10)  three  different  types  of  subdomains: 

Dj  =  conductors,  carrying  currents  impressed  by  external  or  internal 
voltage  or  current  sources  (conductors  stranded  and  transposed); 

D£  =  solid  conductors  which  are  seats  of  eddy  currents; 

=  nonconductive  domains. 
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Figure  10.  Structure  of  the  electromagnetic  system  analyzed  by  the  3-D, 
transient  f.e.m.  code 


The  general  assumptions  are: 

a.  The  magnetoquasistatic  approximation  applies  to  all 
electromagnetic  launcher  problems: 


b.  The  materials  are  isotropic, 

c  Hysteresis  effects  are  neglected. 

3.3  MAXWELL’S  EQUATIONS 

Maxwell's  equations,  as  applied  to  electromagnetic  launchers,  are: 


43 


divB  =  0 
divl  =  0 

fl  =/i77 

~J  =  a(E  +  Elr*  =  7,  +7J  +  7/+.. 

where  the  "impressed"  fields: 

Eimp(.l  —  1,2,... /t) 

correspond  to  "impressed"  electromotive  forces  of  nonelectric  origin  (due  to 
chemical  effects,  thermal  effects,  etc.). 

In  the  case  of  electromagnetic  launchers,  persistent  currents  induced  in 
the  projectile  in  the  earlier  stages  of  launching,  as  well  as  the  electromotive 
forces  of  electromechanical  origin  obtained  by  integrating  the  (uxfi) 
induction  terms,  must  be  taken  into  account  properly  as  being  "impressed." 

3.4  GOVERNING  EQUATIONS 

A  theorem  due  to  Helmholtz  [25],  [32]  states  that  any  vector  field  which 
is  finite,  uniform,  and  continuous  and  which  vanishes  at  infinity  may  be 
expressed  as  the  sum  of  a  gradient  of  a  scalar  V  and  a  curl  of  a  zero  divergence 
vector  A- 

F  =  grady f  +  curlA 


(3.4) 

(3.5) 

(3.6) 

(3.7) 


under  the  condition: 
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divA  =  0  (3.9) 

then,  this  separation  into  a  divergenceless  part,  curl  A ,  and  a  curless  part,  grad 
y '  is  unique  [25],  [32]: 

B  =  curlA  (3.10) 


The  Coulomb's  gauge  applies  in  order  that  the  continuity  equation  be 
satisfied  in  the  domain  of  interest.  It  is  enforced  in  the  f.e.m.  formulation  by 
use  of  Lagrange  multipliers,  as  shown  in  the  next  chapter. 


Faraday's  law  is  expressed  (3.3)  as: 


V  x  E  = 


dB 

dt 


However,  when,  in  a  body  (conductor  or  not),  the  domain  in  which  the 
field  £  is  calculated  moves  with  velocity  “  (as  in  the  case  of  the  projectile  of 
the  electromagnetic  launcher),  then  the  derivative  of  the  magnetic  flux  must 
be  considered  in  a  context  of  an  integration  surface  moving  with  the  material 
[25]. 


Applying  the  general  formula  of  the  flux  denvative  [25]  to  the 
Faraday’s  law  (3.3),  we  obtain: 


—  +  udivB  -  curllu  x  B ) 
dt  x  ' 


•  ds 


(3.11) 


and,  since  divB  =  0,  we  obtain: 
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-  dB  , 
cur  IE  =  — —  +  curllu  x  B 
dt  v  ' 


and,  using  (1.9) 

curlE  =  curlA  +  curl(u  x  Bj 


(3.11a) 


(3.11b) 


Since  the  inversion  of  the  order  of  the  space  and  time  derivatives  is 
permissible. 


cur  IE  =  -cur/j^A  +u  xfij 


(3.11c) 


curl^E  +A+uxB  j  = 


(3.11d) 


meaning  that  the  expression  in  the  parentheses  is  curlless  and,  consequently, 
may  be  expressed  as  a  gradient  of  a  scalar  potential 


E  +  A  +  u  x  B  =  -gradys 


(3.12) 


We  may  consider  a  "generalized"  time  derivative  [25]  of  the  vector 
potential: 


A =  A+uxB 


(3.12a) 


such  that 


E  +  A  =  -gradyf 


(3.12b) 


In  fact,  grady  may  contain  different  te^ms 


46 


grad  y  -  gradyx  +  grady2+... 


(3.12c) 


Then,  the  current  density  is: 

J  =  alE  +  E^)  =  o^-A-grjdy/)  +  ls  =  Jt+7s 


(3.13) 


By  using  (3.13),  the  magnetic  circuit  law  (3.2)  may  be  written  as: 
curlH  =  7  =  o{-  A  -  grady^  +  J 


(3.13a) 


and  by  rearranging  the  terms,  we  obtain: 
<7  A  +  cgrad  y  +  curl — cur l A  =  Jt 


(1.13b) 


In  the  stator,  and  in  the  projectile,  as  well  as  a  corollary  of  Coulomb's  gauge, 
the  current  density  is  solenoidal,  (3.5). 


Thus,  by  considering  equation  (3.13),  we  can  write: 


div  | 


dA  j  ^  7 
--grady  +/, 
dt 


=  0 


(3.14) 


Using,  from  now  on  as  shorthand,  the  del  (V)  operator  and  considering,  also, 
that  A  implies  A ,  we  can  write: 


V* 


(  dA  ^ 

-a — 


+  V*(-<tV^)-V»7j  =  0 


(3.14b) 


Summarizing  the  mathematical  transformations  performed  above,  the 
finite  element,  three  dimensional  formulation  of  the  transient 
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electromagnetic  problem  is  governed  by  the  following  system  (I)of  differential 
equations  with  partial  derivatives: 


a +  V  x  —  V  x  A  +  crV^-  =  J 
dt  (j. 

(3.13b) 

(  d~A\ 

V«  -a —  +  V»(-ctVi/a)  =  0 

V  dt  ) 

(3.14b) 

V*  A  =  0 

(3.9) 

v*j  =  o 

(3.5) 

In  electromagnetic  accelerators  and  their  power  supplies  stranded  and 
transposed  conductors  are  widely  used.  Such  a  situation  is  characterized  by 
the  suppression  of  the  eddy  currents  which  would  otherwise  perturb  the 
quasistationary  distribution  of  the  current  densities  in  the  conductors.  In  the 
system  of  equations  (I),  such  conductors  (shown  as  subdomains  &j,  >&jm) 

are  considered  as  having  zero  conductivity,  since  imposing  cr  =  0  in  the 
numerical  code  forces  the  eddy  currents  J ,  to  be  zero.  Then, 


is  prescribed  as  a  function  of  time,  for  example,  and  the  condition  (1.14d) 
becomes: 


v./;  =0 


(3.5a) 


Using  circuit  analysis  methods,  J ’  can  be  obtained  with  good  approximation. 
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3.5  THE  FINITE  ELEMENT  FORMULATION  FOR  THE  3-D,  TRANSIENT 

ELECTROMAGNETIC  PROBLEM  USING  GALERKIN  METHOD 

The  Galerkin  method  [33],  [34],  [10]  permits  one  to  obtain  an 
approximate  solution  to  a  differential  equation.  In  the  finite  element 
method,  the  application  of  Galerkin  technique  is  based  on  the  requirement  of 
orthogonality  between  the  functions  used  for  approximation  and  the  error 
between  the  approximate  and  the  exact  solution  to  the  differential  equation. 
Considering  a  differential  equation  Lu-  f  ~  0  (where  L  is  a  linear  differential 
operator)  which  has  an  approximate  solution  u  then  the  error  (or 

residual)  e  =  Lu-  f  can  be  efficiently  minimized  if  we  require  the  annulment 

of  integrals  jN;£dw  =  0  for  each  0f  t^e  basis  functions  Nj. 
o 

These  sets  of  integrals  state  mathematically  that  the  basis  functions 
must  be  orthogonal  to  the  error  over  the  entire  domain  D. 

The  finite  element  procedure  is  deeply  rooted  in  variational  calculus, 
the  solutions  being  obtained  by  minimizing  a  functional  over  the  domain 
under  consideration.  The  Galerkin  method  permits  that  finite  element 
formulations  to  be  extended  to  domains  where  no  variational  principle  is 
feasible.  The  transient  electromagnetic  problem  of  diffusion  of  fields  is  one  of 
such  problems  which  being  described  by  parabolic  differential  equations  with 
partial  derivatives  refer,  to  lossy  systems. 

The  finite  element  formulation  for  the  3-D  problem  looks  at  first  sight, 
similar  to  the  formulation  used  for  the  much  simpler  problem  of  two- 
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dimensional  axisymmetric  diffusion  (Pillsbury  [24]).  The  simplicity  of  the  last 
one  is  assured  by  the  need  to  find  only  one  component  for  the  vector 
potential  [a  =  Ae]  while  in  the  three-dimensional  problem  which  is  the  object 

of  the  present  work  it  is  necessary  to  find  all  three  components  of  the  vector 
potential  [a  =  A1  +  Ay+ A,  J . 

As  in  the  simpler  problem  to  the  unknowns  (A  =vector  potential, 

¥  =scalar  potential,  and  ^  =Lagrange  multiplier),  we  associate  a  set  of  test 
functions  u,  v,  and  r,  respectively. 

Among  the  conditions  required  in  applying  the  Galerkin  method  are: 

a.  The  test  functions  must  be  chosen  from  the  same  family  as  the 
tried  functions; 

b.  The  test  and  the  trial  functions  must  be  linearly  independent; 

c  The  trail  functions  should  satisfy  the  boundary  conditions  and 
initial  conditions  exactly. 

Condition  a)  defines  the  Galerkin  method,  condition  b)  is  necessasry  to 
ensure  that  independent  equations  are  available  to  obtain  the  unknown 
coefficients;  and,  finally,  condition  c)  relates  to  the  efficiency  of  the  method. 

In  order  to  apply  the  Galerkin  method,  the  dot  product  of  the  test 
functions  with  the  governing  equations  is  integrated  over  the  volume  ( w )  of 
each  element  and  summated  over  all  k  elements. 


~r 


f  a A+Vx  — VxA  +  crVy/--/. 


<udw+  J/l(V«  u)ftv  +  J  r(V»A)dw 

Dt  Dk 


^  =  0 
Vu,-r 


jfv*^crAj  +  V«(-crVi//-)jwivv 

d,l  '  1 


=  0 


Vv 


Using  a  vector  identity,  the  element  equation  is  transformed,  yielding 
expressions  which  contain  also  the  contributions  over  the  boundary  of  the 
elements  SDk. 


V  x  A 


\ 


1.0  ‘ 


ff— V  x  A 

DJ>  J 


udw  =  f  V  •(  —V  x  A  I  x  udw  +  [ 

i  U  J  i 

(V  x  u)dw  +  f  j  —V  x  A  X  i?  •ruft  = 

SDt  ' 

-  (V  x  u)dw  +  J (//  x  u)  •  lids  - 


—V  x  A 


(V  x  u)  pw  = 


iD, 


ff— VxA  |•(Vxi^)dw’+  |  («  x  H)»uds  = 

D.V^  /  SZ3, 


where  the  integrals  over  Djt,  are  volume  integrals  and  over  SDk,  ar^  surface 
integrals,  respectively. 


In  the  same  manner: 
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J^crA  j•(Vv)dw-  Jv(-ctA  j*  nds  = 

Dk  SDk 

and  also: 

J[7  •  (-cr  V  \f/)^dw  = 

Dk 

J(cr  Vi^)(Vv)dw+  Jv(-ctVi j/)nds 

Dk  SDk 

If  the  integrals  referring  to  the  test  function  v  are  summarized,  their 
physical  meaning  becomes  evident: 

j  jv  •  (— cr  a)  +  V  •  (-or  V^)jvdw  = 
d}  i 

jf|crAj»(Vv)  +  (o'V^)«(Vv)W+  Jv^-crA  - aV^jnfds  = 

Dk  ■*  SDt 

|  o|a  +  V  i^jVvdw  +  |  vJt  •  rids  = 

Dk  SD„ 

J-74Vvdw  +  J  Jt  •  nvds  =  0 

Dk  SDk 

Vv 

Finally: 

J  cr  A*u  +  — VxA  «(V  xu)  +  a  Vy/-*  u -J,  •  «  +  A(V«u)+  t(V*  a)  dw  + 
okl  U  ) 

J(n  x  H)»uds  =  0 

SDk 


Vm,t 
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and 


JaU+Vvr]«Vvdw+  J Jt 

Dk  SDt 


•  nvds  =  0 


Vv. 


The  trial  functions  representing  the  magnetic  vector  potential  a  ,  the 
electrokinetic  scalar  potential  V ,  the  Lagrange  multipliers  A  and  the 
corresponding  test  functions  u,v,  and  t  must  be  smooth  enough  in  order  to 
assure  that  the  integrand  in  the  Galerkin  formulation  is  finite.  If: 

A,y,urv&  H1 

and 


A, re  L2 


Then  the  weak  form  obtained  by  using  the  Galerkin  technique  leads  to  the 
finite  element  equations  written  on  element  basis: 


i 


a  A  *u  -—V  xA»Vxu  +  dV\jfu-Jsu  +  A(V»u)+  r(V  •  a) 

r* 


21 

*  o, 

^  J  [i^x  Hl»H[+'n2H2»n2\cls+  jnxH  •nds  =  0 

*  SDk-SD  SD 


aw  + 


Vu,T 

For  electromagnetic  accelerators  and  their  power  supplies,  a  minimum 
electromagnetic  signature  represents  a  figure  of  merit. 


53 


Then,  prescribing  a  zero  value  for  the  potential  vector  at  far-field 
( R  =  °°)  represents  a  normal  condition  and  occurrence.  The  values  of  the 
fields  become  negligibly  small  at  relatively  short  distances  because  of 
compensations  and  efficiency  considerations.  It  follows  that  at  the  boundary 
SD  of  the  domain  being  analyzed,  u  =  0,  and  this  means  that: 

J(n  x  77)  •  uds  -  0, 

SD 

The  assumption  was  made  that  u  e  implying  that  u  is  continuous 
across  the  interface,  then: 

^  Jf/i,  x//,  •iij  ■= 

*  SDk  -SD 

X  J["i  x(^i~W2)]«  uds  =  0 

*  SDk-SD 


since 


u,  x  //,  =  n,  x  H2 


As  a  result,  the  finite  element  equations,  written  on  the  element  basis 
using  a  Galerkin  technique,  simplifies  to: 


X 


o. 


crA»u  +  -VxA*Vxu  +  (jVw«il-7*u 


UW  + 


y=o 


J[-A(V*u)+  -r(V*/l)|div 


(3.16) 


The  "electrokinetic"  problem  considerations  lead  to  the  finite  element 
equations  written  on  element  basis: 
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^  J  c|a  —  V i/aJ  •  Vvdw  +  ]|T  +  Jti  •  n^Jfls  +  J/,  •  nvds  =  0 

*  Dt  SDt  -SD  SD 


Vv 


Using  the  same  far-field  considerations  as  above,  we  assume  with 
sufficient  accuracy  that  the  scalar  (electrokinetic)  potential  is  zero  at  the 
boundaries  of  the  domain  (v=0  on  SD).  Then, 

J  Jt  •  nvds  =  0 

SD 

By  virtue  of  the  assumptions  made  before,  u  e  H1,  v  is  continuous  at 
the  interfaces;  then: 

SDt-SD 

SDk-SD 

since  7 ,  •  n,  =  7^. 


then: 


J|crA  •  Vv  +  a  V y/ •  Vv^iw  =  0 

k 


(3.17) 


Relations  (3.16)  and  (3.17)  represent  the  final  form  of  the  system  of 
equations  on  element  basis  for  the  Galerkin  approach  of  the  3-D  transient 
finite  element  code. 
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In  actually  solving  a  particular  electromagnetic  problem,  it  is  important 

to  systematically  look  at  the  form  taken  by  the  boundary  conditions  for  H,  B, 

and  /: 

n  x  H  =0; 

7  •  n  =  0; 

B  •  n  =  0; 
n  x  J  =  0; 

For  the  potentials  (vector  and  scalar)  A  and  yr/  a  similar  systematic  look 
is  necessary: 

V«(nxA)  =  A»Vxn-n«VxA  =  n*B 
If  n  x  A  =  0,  then  -B  •  n  =  0 

nx/=nx  |  <7  A  -  <rVy/j  = 

-o|n  x A  +  nxVi/Aj  =  «xVy/ 

For  nonconductive  materials,  as  well  as  for  stranded  and  transposed 
conductors,  <r  is  considered  zero: 


<7  =  0  nxj=  0 

For  other  conductive  materials  cr  *  0,  but  since 

V  ~  V  ~  constant,  then  nxj  =  0  (where  i  =  2,  2..n)  n  being  the  number  of 


conductors. 
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3.6  FINITE  ELEMENT  INTERPOLATION 
I  =  Afo 


u  =  uiei 

y/  =  Nly/t  A  =  M"Aa 


v  =  ?=MmTM 

dA  dA  dAe  dAfi 

*  CT_.  w  =  rr-~Luj8t/  =  cStjNp  ~Nau“  =  oN^pd^uf 


*  crVy/  •  u  =  ay/  ^  =  crA*  =  crA^Na 


*  y  •  w  =  y^x 

*  —V x  A  •  VxET 


H  ep<rAr,qepsiUt,s 

=  ^(5Jrt-SqA,)Ar.qUl.1 
=  ^(A'.sU<,s  ~  As.'U,.s) 

= j(n,,Xn«X  -  n,Xn°X) 
=  \(ns,N„X<  - 

r*’ 


where  erst  is  the  permutation  symbol  defined  by  the  equation: 

ein~e222=e333=e112=e113=  "=0 

e123=e231=e321=1 

e213=e321=e132='1 
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equal  to  zero  whenever  the  values  of  any  two  indices  coincide;  eikj  equals  1 
when  the  subscripts  permute  like  1,2,3,  and  Cyjt=-1  otherwise. 

If  8  is  the  Kronecker  symbol,  e  relates  to  <5  according  to  the  identity: 
=  sjAsj'S* 

*  A(V  •u)  =  X  Uiit  =  M,  XHNaX  = 

*  t(V  •  A)  =  TAJtJ  =  MmNpjTM  A1} 


3.7  ELEMENT  EQUATIONS 


Vvt 


Then: 
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J«rW„V* 


A 


dA f 
”dT 


■+ 


.  J;(WrV») 


J  aVX 

A 


¥t  + 


j  M"N„„ 


where  i  =  1,2,3  and  a  =  1,2 ,...N,  N  being  the  number  of  nodes. 

dA„  dA.  dA?  R  k dA f 

*  cr — Vv  =  <r — -v  =  <jiVB — J-Nfivk  =  (tNbN, — -vt 

dr  dr  '7  pdt  p  '*  dt 


*  <rVy/Vv  =  ctyiSvs  =  oN;  ytN\vk  =  a N‘tN,  Yivt 


MmN , 


A?  =0 


m  =  1,2,...  M 


where  M  is  the  number  of  multipliers 


j W, 

0, 


<£4? 

dr 


j  on'n. 


Ld‘ 


V/  =0 


We  denote: 


X  = 


4 

Vt 

X. 


x  = 


dAp 


dt 

dyt 


dt 

dX. 


dt 


(3.18) 


Then,  the  system  of  equations  in  compact  matrix  form  is: 
Cx+kx  =  F 


(3.19) 


where 
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Pj  t  n 

\aNaNp8ljf  0,  0 

«  Dk 

C  =  k  ](TNPN^  0,  0 

m  0,  0,  0 


(3.20) 


(3.20a) 


3.8  TIME  MARCHING  ALGORITHM 


(3.20b) 


In  solving  the  transient  (diffusion)  electromagnetic  problem  described 
by  the  finite  element  system  of  equations: 

Cx+  kx  =  F 

m  -  »  —  - 


(3.19) 
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the  first  derivative  of  the  nodal  point  unknowns  is  approximated  by  the 
discrete  form: 

jc  =  (1  -d)xn+ex^x 


.  xn+1-xn 

x  = - 

At 

leading  to: 

+  A:[(l-0)x"+6x"+1 

\  J 


=  (1  -Q)Fn+dF 


*+ 1 


and  further: 


'  c 

r 

-=-  +  6k 

x"+1  = 

At 

At 


-{l-6)k 


jc"+(1  ~e)Fn+dF 


when  0  =  1  we  have  a  backward  difference  form  of  the  ordinary  differential 
equations  which  is  unconditionally  stable. 


The  matrix  equation  becomes: 


CA 

I  1-A  lc¥  kx 

'  "'OBJ 

^  0 
e  kx  o,  o 


$ 

¥t 

A 


'(«») 


C^',  0,  0 

A* 

f  J*"Na 

A  t 

i 

o, 

a 

> 

o 

> 

o 

• 

Vt 

0 

0,  0,  0 

K 

0 

(3.21) 
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where 

CA  = 

l<rNJi,S, 

£ 

II 

Dk 

\cNtN'J 

II 

ok 

j-jKAA 

Dkr 

II 

-i* 

\  °w. 

lcx  = 

^CUn 

Dk 

|  M"N„, 

Sf^ 

II 

Dk 

J  °KK 

II 

Dk 

\m-ns, 

Dk 

o 


(3.21a) 
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CHAPTER  4 

SUPERFICIAL  CURRENTS  AS  LAGRANGE  MULTIPLIERS  IN  THREE- 
DIMENSIONAL  FINITE  ELEMENT  METHOD  FORMULATIONS  FOR 
ELECTROMAGNETIC  PROBLEMS 

4.1  INTRODUCTION:  THE  ENFORCEMENT  OF  COULOMB'S  GAUGE 

f— =ol 

In  conditions  of  quasistatic  approximation  ^  dt  J  which  apply  to  all 

electromechanical  magnetic  flux  compression  and  magnetic  flux  expansion 
converters.  Coulomb's  gauge  [31],  [32],  [10]  for  the  magnetic  vector  potential  is 
"chosen:" 

divA  =0.  (4.1) 

The  gauge  "choice"  is  necessary,  since  a  vector  field  must  be  completely 
characterized  by  specifying  both  its  curl  and  its  divergence  in  every  point  of  its 
domain  of  definition  [25],  [31],  [32].  While  drastically  simplifying  the 
governing  equations  (see  Chapter  3),  the  main  reason  for  introducing  is  not  to 
simplify  them,  as  stated  in  [24],  but  is  enforced  in  order  to  satisfy  the 
continuity  equation  in  the  domain  of  interest: 

divJ  =  0.  (4.2) 

Equation  (4.2)  is  actually  a  corollary  of  the  Coulomb  gauge  for  the 
condition  that: 


r2\a\ <  ~ 


when  R 


(4.3) 
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since,  by  taking  the  divergence  of  the  two  members  of  the  Helmholtz 
equation  for  magnetic  fields  (simplified  for  the  magnetostatic  conditions  [25], 
[31]): 

V2J  =  -fi  •  J ,  (4.4) 

and  considering  (4.2),  we  find  that  divA  is  a  constant  and,  since  its  value  is 
zero  at  the  infinity,  the  constant  is  zero  everywhere. 

The  enforcement  of  the  Coulomb  gauge  will  mean  that,  in 
electromechanical  converters,  all  the  conduction  currents  are  closed  in 
themselves,  and  the  energy  spent  in  electromagnetic  radiation  is  negligibly 
small,  even  for  the  case  of  rapid  launchers. 

For  electromagnetic  systems  modelling,  there  exist  several  two- 
dimensional  f.e.m.  numerical  codes.  Their  usefulness  is  based  on  the  fact  that 
many  electromagnetic  problems  can  be  simplified  to  two-dimensional 
symmetries  and  solved  with  more  modest  computational  means. 

In  such  codes,  the  equation  (4.2)  and  the  Coulomb  gauge  are  easily 
enforced.  In  an  axisymmetric  structure,  a  coil  appears  as  its  cross-section  and 
is  implied  that  the  current  distribution  will  be  identical  in  ail  its  other  cross- 
sections. 

In  three-dimensional  structures  modelling  the  true,  unsimplified 
structure  of  the  electromagnetic  system,  the  situation  changes,  especially  in 
the  finite  element  method  since  the  entire  domain  is  built  by  volume 
elements  (bricks,  pyramids,  wedges,  etc.)  which  must  match  the  continuity 
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equation  (4.2)  not  only  in  some  points  of  the  integration,  but  along  and  all 
over  the  entire  domain.  In  other  words,  the  special  role  of  enforcing  the 
Coulomb  gauge  in  three-dimensional  f.e.m.  systems  is  to  bridge  a  method 
which  is  fundamentally  discrete  with  a  requirement  of  pointwise  continuity. 
The  "plaster"  used  to  "smooth-out"  the  discrete  3-D  electromagnetic  structure 
is  a  group  of  superficial  currents  flowing  through  the  3-D  element  faces,  thus 
producing  magnetic  fields  which  impose  the  constraint  divA  =0.  Such 
superficial  currents  will  be  assimilated  with  Lagrange  multipliers. 

This  chapter  will  be  devoted  to  outlining  the  special  physical  meaning 
of  superficial  Lagrange  multipliers  in  three  dimensional,  finite  element 
formulation  of  electromagnetic  problems,  in  magnetic  flux  compression 
devices. 

4.2  THREE-DIMENSIONAL  ELEMENTS 

Finite  element  method  [22],  [28],  [29],  [30]  is  based  on  the  replacement  of 
a  continuous  system  (which  ,  in  our  case,  is  the  domain  in  which  the 
electromagnetic  transient  problem  must  be  solved)  by  an  equivalent  discrete 
system,  made  out  of  finite  elements  (volume,  or  three-dimensional,  in  our 
case).  The  geometry  of  all  the  elements  :  'ust  be  defined  analytically  and 
appropriate  interpolation  functions  must  be  constructed  for  each  element. 

To  define  the  geometry  of  elements,  a  set  of  n  points  is  selected.  These 
points,  called  geometrical  nodes,  may  sometimes  coincide  with  interpolation 
nodes.  Each  element  is  analytically  and  uniquely  defined  in  terms  of 
geometrical  nodes  belonging  to  that  element  and  its  boundary.  Two  distinct 
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elements  can  have  common  points  only  on  their  common  boundaries  if  such 
boundaries  exist;  no  overlapping  is  allowed.  Figure  11  shows  elements  as 
rectangular  prisms  (linear  and  quadratic)  with  8  and,  respectively,  20  nodes. 

In  order  to  assure  the  interelement  continuity,  it  is  necessary  that,  along  a 
whole  face  of  an  element,  the  nodal  values  define  a  unique  variation  of  the 
unknown  function. 

Using  three  normalized  coordinates,  we  have  the  following  shape 
functions  [27],  [30]: 

For  the  linear  element  (8  nodes): 

=i(1+<’°X1+ 7i°)(1+C°) 

V  (4.5) 

For  the  quadratic  element  (20  nodes): 

Comer  nodes: 

n,  =  i(i+ ?o)(l + »fe)(l + f.X$>  +  no  +  Co  -  2) 

8  (4.6) 

Mid-size  nodes: 


|,=on,=±if,=±i 

(47) 

Figure  12  shows  triangular  prism  elements,  linear  and  quadratic,  with  6 
and,  respectively,  15  nodes.  In  this  case,  the  shape  functions  are  obvious  and 
found  in  any  book  on  finite  element  method  [27],  [30].  Again,  for 
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Figure  11.  Three-dimensional,  finite  elements  as  rectangular  prisms  with  8  or 
20  nodes.  Surface  currents  through  faces  represent  Lagrange  multipliers. 
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T|  =  1 


Figure  12.  Three-dimensional,  finite  elements  as  triangular  prism  elements 
with  6  or  15  nodes.  Surface  currents  through  faces  represent  Lagrange 
multipliers. 
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interelement  continuity,  it  is  necessary  that,  along  a  whole  face  of  an  element 
the  nodal  values  define  a  unique  variation  of  the  unknown  function. 

4.3  PRESCRIPTION  OF  THE  FORCING  FUNCTION  IN  THREE- 

DIMENSIONAL  FINITE  ELEMENT  ELECTROMAGNETIC  PROBLEMS 

In  electromagnetic  problems,  the  forcing  functions  are  globally,  voltage 
sources  or  current  sources,  impressed  upon  a  certain  network.  In 
electromagnetic  field  problems  solved  by  f.e.m.  codes,  the  forcing  functions 
corresponding  to  voltages  are  the  intensities  of  applied  electric  fields  and, 
corresponding  to  current  sources  are  current  densities  . 

Current  densities,  for  example,  are  prescribed  element  by  element. 

Two  distinct  elements  have  as  common  boundaries  points,  lines,  or  surfaces. 
At  their  common  boundaries,  the  assembled  elements  should  leave  no  holes 
or,  as  mentioned  above,  overlap.  When  the  boundary  of  the  domain  cannot 
be  exactly  filled  by  the  elements  selected,  an  error  called  geometrical 
discretization  error,  cannot  be  avoided. 

By  prescribing  current  densities,  element  by  element,  the  current, 
representing  the  surface  integral  (the  flux)  of  the  local  current  densities,  does 
not  obey  exactly  the  first  Kirchhoff's  law  or  locally  divJ  =  0  (or  by  implication 
divA  =  0)  does  not  hold,  two  elements,  one  hexahedronal,  the  other  prismatic, 
may  be  adjacent  and  the  current  density  prescription  as  forcing  function  may 
lead  to  discrepancies  in  the  Coulomb's  gauge  which  must  be  fulfilled 
pointwise. 
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Round-off  errors  in  computation  of  currents  -  even  for  an  apparently 
balanced  input  -  lead  to  discrepancies  upsetting  the  solenoidal  character  of  the 
current  densities. 

The  method  of  correcting  such  discrepancies  is  the  introduction  of 
surface  (sheets)  currents  flowing  through  the  element  faces  and  by  the 
magnetic  field  they  create  a  compensating  effect.  Then,  the  Coulomb's  gauge 
seems  to  be  imposed  pointwise  in  the  entire  volume  of  the  domain.  We  are 
identifying  such  sheets  of  currents  with  Lagrange  multipliers. 

4.4.  PHYSICAL  INTERPRETATION  OF  SUPERFICIAL  CURRENTS  AS 

LAGRANGE  MULTIPLIERS  IN  F.E.M.  FORMULATIONS  OF  THREE- 

DIMENSIONAL  ELECTROMAGNETIC  PROBLEMS 

The  Helmholtz  theorem  [25],  [31]  states  that  any  vector  field  is  uniquely 
separable  into  a  divergenceless  part,  curl  A  and  a  curless  part,  grad  W  t  in  each 

point  of  the  field,  and  in  average  in  each  of  the  three  dimensional  elements 
in  the  finite  element  formulation.  We  are  interested  in  the  transverse 
component  of  the  vector  field  (transverse  to  the  direction  of  greatest  change 
which  is  the  lamellar  or  longitudinal  one),  the  Coulomb  gauge  imposing  a 
linear  relationship  between  the  three  components  of  the  vector  potential  A , 
decreasing  the  number  of  independent  components  to  two. 

For  the  3-D  f.e.m.  formulation  of  the  transient  electromagnetic 
problem,  the  Galerkin  method  is  used  in  order  to  obtain  the  finite  element 
equations.  The  formulation  is  carried  out  in  Chapter  3.  The  equation  (3.18) 
from  this  chapter  is  reproduced  here: 
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i  cr— +  Vx  — xVx  A  +  crVu/--/. 
U  dt  V. 


•  u  + 


A  ( V  •  u )  +  t(V  •  A  )}dv  =  0 


(3.18) 


for  any  u  and  x- 

In  equation  (3.18),  the  Coulomb  gauge  constraint  div  A  =0  has  been 
enforced  using  the  Lagrange  multipliers  A  and  the  corresponding  test 
functions  t. 

From  the  above  equations,  the  dimensionality  of  the  Lagrange 
multiplier  results  in  and  represents  a  surface  current  (a  current  sheet). 

In  any  type  of  3-D  elements  (as  in  Figure  11,  for  example),  they 
materialize  in  the  face  of  the  elements.  Their  use  is  based  on  the  fact  that  the 
magnetic  field  produced  by  a  distribution  of  volume  currents  inside  an 
element  (brick,  prism,  tetrahedron,  etc.)  can  be  substituted  by  an  identical 
magnetic  field  structure  produced  by  a  distribution  of  surface  currents  flowing 
on  different  faces  of  the  respective  element. 

The  current  density  for  the  sheet  current,  K  is  obtained  by  passing  to 
limit  of  a  volume  current  density  J  when  one  of  the  dimensions  (the 
thickness  t )  goes  to  zero  (Figure  13). 

The  total  current  I  (Figure  14)  can  be  expressed  as 
/  =  7  •  n(lt)  =  K  •  nt(t) 


Figure  13.  Surf; 
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Figure  14.  Cal 
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where 


K  =  limi  •  t 


(4.8) 


In  general,  the  surface  current  7  (Figurel4)  is  given  by  the  integral 


N 

I  =  \(K»nt)dl 

M 


(4.9) 


The  surface  currents  equivalent  of  the  continuity  condition  for  the 
volume  currents  expressed  in  point  form  by  equation  (4.2)  is: 


j(K*nt)d£  =  0 


(4.10) 


In  some  treatments  (32],  equation  (4.8)  is  replaced  by  a  double  integral  of 
a  surface  divergence  of  the  current  density  K . 

To  illustrate  how  the  LagTange  multipliers  act  in  imposing  a  constraint 
(e.g.  the  gauge  condition),  let's  assume  a  given  distribution  of  currents 
J(x,y,z)  in  a  volume  V  already  divided  in  3-D  elements.  The  volume  V  is 
surrounded  by  the  closed  surface  £  such  that  the  flux  density  B,  produced  by 
the  volume  currents  J  is  tangential  everywhere  to  the  surface  £ . 
Adistribution  of  sheet  currents  of  density  K  over  the  surface  Z,  is  chosen 
(Durand  [32])  such  that: 


-IL0K  =  n  x  B, 


(4.11) 
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Then,  multiplying  by  the  positive  (outward)  normal  to  £,  both  sides  of 
(4.11)  and  considering  the  properties  of  the  double  cross  product,  yields  the 
value  of  the  flux  density  results  as: 

B,=fX0(nxK)  (4.12) 


If  we  superimpose  (considering  linear  media)  the  effects  of  the  two 
current  distributions,  volume  J  and  surface  K  we  find  the  flux  densities: 


B=B, 

inside  of  surface  I 

B=0 

outside  of  E . 

The  known  condition  of  discontinuity  (jump)  in  the  value  of  B  when 
crossing  a  sheet  of  currents 

B above  -  B  below  =  /*«["  x is  satisfied  by  the  relations  (4.12)  and  (4.13). 


From  Biot  and  Savart’s  law: 

* = H(*  x  ?)* = rt* x  grad^)Y = 


(4.14) 


and  knowing  both  its  divergence  and  its  curl  ( divB  -  0,curlB  -  ^i0J),  B  is 
completely  defined  for  the  electromagnetic  problems  concerning 
electromechanical  converters  (launchers  and  electromechanical  power 
supplies.) 


74 


If  the  currents  K  had  existed  alone,  they  would  have  produced  outside 
of  surface  £  a  flux  density  equal  and  opposite  to  that  produced  by  the  volume 
distribution  of  currents  Jand  a  flux  density  B  equal  to  zero  inside  of  £. 

On  this  fact  lies  the  physical  meaning  and  the  mechanisms  of 
operation  of  the  Lagrange  multipliers,  acting  as  current  sheets  in  the  3-D 
elements  in  a  finite  element  grid:  They  are  the  means  to  replace  the  currents 
7  for  their  action  outside  the  boundary  £  by  superficial  currents  flowing  over 
surface 

In  the  same  way,  we  can  consider  the  case  in  which  the  closed  surface  £ 
is  exterior  to  the  volume  enclosed  by  the  3-D  element  (or  the  set  of  3-D 
elements  of  which  the  volume  of  definition  for  the  electromagnetic  problem 
consists).  Following  the  same  reasoning  as  above,  the  superficial  distribution 
of  the  sheet  of  currents  K  produces  inside  of  boundary  £  the  same  flux 
density  °it 

B,  =  -il\n  x  K\  (4.15) 

as  the  volume  currents  7.  Outside  of  surface  2,  the  flux  density  produced  by 
the  current  sheet  K  is  zero.  As  an  example,  in  Figure  15  we  have  a  cylinder 

with  a  radius  r  through  which  a  current  I  of  density/  =  — flows  uniformly. 

Jtr 

On  the  toroidal  surface  surrounding  concentrically  the  cylinder,  a 
surface  current  density  K  can  always  be  chosen  in  order  to  produce  a  field 
inside  of  the  torus,  identically  equal  to  zero. 


Figure  15.  Illustrative  of  the  mechanism  of  current  sheets  acting  as  Lagrange 
multipliers. 
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We  must  emphasize  that  Lagrange  multipliers  refer  to  the  vector 
potential  A,  enforcing  its  divergence  to  be  zero,  not  to  the  flux  density  B 
which  is  the  curl  of  A  -  However,  since  the  topic  is  the  physical  significance  of 
the  Lagrange  multipliers,  the  flux  density  which  is  more  familiar  was  used. 


Lagrange  multipliers  impose  the  condition  div  A=0.  The  expression  of 
divergence  in  terms  of  surface  currents  for  the  magnetostatic  problem  results 
from  [25],  [32]: 


d‘VA  = 


(4.16) 


and  since  K  is  a  vector  tangent  to  the  surface,  we  may  write  equation  (4.16)  as: 


divA =fe{ndiv-*7"Jidiv'(7)* 


(4.17) 


Where  the  following  vectorial  identity  was  employed: 


div  —  ]  =  -divsK  -  Kgradi- 
r  r  \  \r 


(4.18) 


In  (4.17)  and  (4.18),  the  s  index  was  used  in  order  to  indicate  the  surface 
gradient  and  divergence. 


Since 


jjdiv,Kds  =  \(n*K}it  =  0 


(4.19) 


where  the  curve  T  surrounds  the  surface  2,  we  conclude  that: 
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divsK  -  0 

Also,  on  the  curve  r  we  have: 

0 

resulting  that  =  0  is  imposed  by  surface  currents  flowing  on  the  faces  of 
elements  in  electromechanical  converters  treated  here. 

The  boundary  conditions  for  the  vector  potential,  when  crossing  a 
sheet  of  currents  are: 

A- a  ~  \ 

and 

( (dA^  - 
-j-  "  T"  =-VoK 

\dn  Ja  \dn  Jb  (4.20) 

4.5  AN  ILLUSTRATATION  OF  THE  USE  OF  THREE-DIMENSIONAL 
F.E.M.  TRANSIENT  ELECTROMAGNETIC  CODE 

As  an  illustration  of  the  application  of  the  three-dimensional,  finite 
element  transient  electromagnetic  code  using  superficial  currents  as  Lagrange 
multipliers  -  a  simple  flux  compression  structure  was  chosen.  This 
arrangement  comprises  a  rectangular  coil  of  two  turns  made  of  stranded  and 
transposed  copper  conductor  in  which  a  current  is  impressed  in  a  ramp 
variation.  The  current  density  increases  linearly  in  the  coil  from  0  to  6.2  x 
A 

109~T  in  a  time  of  25  gsec. 
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In  front  of  the  coil,  compressioning  the  magnetic  flux  is  a  rectangular 

1 

piece  of  solid  copper,  having  a  conductivity  of  5.8  x  107  &n,  which  is  the  seat 
of  the  eddy  currents.  Referring  to  Figure  10  in  Chapter  3,  the  coil  is 
representative  of  subdomains  Dj,  the  copper  plate  of  subdomains  De,  and  the 
rest  of  the  space,  of  domains  Do. 

The  structure  is  truly  three-dimensional;  however,  it  can  be  grossly 
approximated  by  a  two-dimensional  axisymmetric  configuration  for  a 
comparison  check. 

Figures  17  and  18  represent  the  division  of  the  entire  domain 
considered  in  the  problem  in  3-D  elements.  The  electromagnetic  (transient) 
diffusion  problem  is  solved  at  time  intervals  of  5  psec  and  plotted  at  intervals 
of  5,  10,  15,  and  20  psec.  The  figures  represent  the  magnitude  of  current 
density  in  its  vectorial  representation  and  the  magnetic  field,  respectively,  in 
upper  face  of  the  copper  plate  and  three-dimensionally  at  5,  10,  15,  and  20  psec 
time. 


It  is  clear  that  division  in  a  relatively  small  number  of  elements  takes 
its  toll  in  the  precision  of  representation.  It  is  also  notable  that  the  accuracy 
remains  high,  the  introduction  of  Lagrange  multipliers  playing  an  important 
role. 
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In  order  to  have  a  term  of  comparison,  a  similar,  but  two- 
dimensional,  axisymmetric  system  was  considered.  The  elements  have 
comparable  size,  in  spite  of  differences  in  shape,  the  coil  is  still  stranded 
and  tranposed,  while  the  massive  copper  plate  has  the  same  conductivity 
as  its  three-dimensional  counterpart.  The  current  density  was  chosen  in 
order  that  the  total  magnetomotive  force  would  be  comparable  for  the  two 
coils.  The  plots  are  still  made  for  5,  10,  15,  and  20  psec  times. 

The  axisymmetric  structure  of  the  simple  flux  compiessor  is  given 
in  Figure  2.  Figures  3-6  show  the  unfolding  in  time  of  the  flux  lines  as  the 
coil  is  pulsed  with  current,  and  eddy  currents  are  developing  in  the  plate. 
This  two-dimensional  case  is  reproduced  from  a  paper  by  M.  D.  Driga  and 
H.  D.  Fair,  "Advanced  Concepts  for  Electromagnetic  Launcher  Power 
Supplies  Incorporating  Magnetic  Flux  Compression,"  [35]. 
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ectorial  representation) 


copper  plate  at  t  =  [0  jXsec  (vectorial  representation) 


copper  plate  at  t  =  1 0  pusec  (3-1)  vectorial  representation) 


3.2-01 


id inal 
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Figure  2h  .  MFC:  current  dens  ity  in  copper  plate  at  1 5  ^see  (vectorial  representation) 
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CHAPTER  5 

AN  APPROACH  OF  CALCULATING  FORCES  IN  ELECTROMECHANIC  AT. 
ACCELERATORS  BY  USING  THE  METHOD  OF  FLUX  DERIVATIVE 

5.1  INTRODUCTION:  POWER  BALANCE  AND  THE  MAGNETIC  FLUX 
DERIVATIVE 

The  operation  of  electromagnetic  macroparticle  accelerators  is  based  on 
the  mechanism  of  conversion  of  electromagnetic  into  mechanical  power1. 
Mechanical  power  is  the  rate  at  which  mechanical  work  is  performed  by  a 
distribution  of  forces  acting  on  the  moving  projectile.  In  the  case  of 
electromagnetic  launchers  such  forces  are  of  electromagnetic  origin.  We  will 
consider  three  ways  of  determining  the  mechanical  forces  of  electromagnetic 
principals  responsible  for  driving  the  electromagnetic  accelerators: 

a.  As  the  distribution  of  body  forces  in  the  projectile  volume  and 
in  the  rest  of  the  accelerator; 

b  As  the  global,  resultant  force  (a  generalized  Lagrangian  force, 
acting  along  a  generalized  coordinate); 

c.  And,  finally,  as  stresses,  using  the  artificial  method  of 

electromagnetic  stress  tensor,  reducing  the  problem  of  electric 
and  magnetic  fields  of  force  to  that  of  an  elastic  continuum. 

1  Characteristic  to  electromagnetic  mac-oparticle  accelerators  is  their  extremely  high  power 
ratings,  for  very  short  time  intervals.  Ihe  energies  involved  are  relatively  low  with  respect  to 
other  high  energy  industrial  processes:  1  kwh=3.6  MJ  When  such  energy  is  delivered  in  2  msec, 
for  example,  it  corresponds  to  a  power  of  1.8xl09w=1.8r:w 
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The  theoretical  tools  which  transform  the  expressions  for  forces  from 
one  form  to  another  are  related  to: 

a.  The  constancy  of  magnetic  flux  as  the  constraint  used  to  assure 
the  proper  transformation  in  order  to  find  the  generalized 
Lagrangian  forces  along  a  generalized  coordinate  (global  or 
integral  characterization);  and 

b.  The  complex  and  subtle  notion  of  flux  derivative  [25]  as  applied 
to  electromagnetic  systems  in  motion.  Again,  using  the 
constancy  of  the  magnetic  flux  as  a  constraint,  we  obtain  all  the 
equivalent  body  forces  of  electromagnetic  origin  (local  or  point 
characterization) . 

Since  the  electromechanical  accelerators  and  their  power  supplies  are 
based  on  magnetic  flux  compression  devices,  it  is  very  significant  that  the 
constancy  of  the  magnetic  flux  constraint  in  the  Lagrangian  formulation,  and 
the  same  in  the  flux  derivative  approach,  is  the  analytical  tool  which  unifies 
the  theoretical  treatment  of  such  devices. 

Macroscopically,  forces  of  electromagnetic  origin  can  be  expressed 
directly  and  elegantly  since  the  electromagnetic  component  of  the  energy  of  a 
thermodynamic  system  can  be  distinctly2  separated  from  other  types  of 
energy.  The  instantaneous  power  balance  for  an  electromagnetic  system  can 
be  written  as: 

2  Electromagnetic  energy  of  a  system  is  that  part  of  the  ’"‘"mal  energy  which  depends  only  on 
the  state  variable*  of  the  electromagnetic  field,  varying  when  they  vary  and  remaining 
constant  when  they  remain  constant. 
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<f>  £(E  x  H)  •  ds  =  $  j  jEJdv  +jj  +  H  •  ^-)i  i> 

so  *  0  (5.1) 

In  certain  conditions  of  media  homogeneity  the  balance  can  be  written 
as: 

f>£(£  x a).  *-JJJsa«.  +  £JJj(£  +  £)i» 

0  D  (5.1a) 

or 

£<f>(E  x  H  )•  ds  =  j  j  j  EJdv  +  jjj\  JJ(zee  +  wm)  dv 

S  V  V 

v 

where 


DE 

eE2 

Dl 

2 

2 

2e 

(5.1b) 

BH 

^2 

B2 

2 

2 

2  u 

(5.1c) 

are  the  electric  and  magnetic  energy  densities. 

Or,  in  other  words,  the  input  of  energy  in  unit  time  in  the  system  (of 
volume  /  and  surrounded  by  the  surface  SB )  expressed  by  the  flux  of 
Poynting's  vector  through  the  area  S„  enclosing  the  system  is  spent  as  joule 

losses  (E  •  J  )  and  also  for  increasing  the  electric  energy  density  (  2  )  and 

(HB  ) 

magnetic  energy  density  \  2  /  °^  *he  system. 
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As  it  is  well  known,  electromagnetic  power  balance  (5.1)  is  obtained 
directly  from  Maxwell's  equations,  after  few  mathematical  manipulations, 
starting  from: 


VxH=/+| 
Vx  E  =  J 


(magnetic  circuit  law) 
(Faraday's  law) 


(5.2) 

(5.3) 


and  by  dot  multiplying  by  E  both  sides  of  (5.1  A)  and  using  the  vectorial 
identity: 


V  •  (e  x  H  )  =  H  •  V  x  E  -  E  *VxH 


we  obtain: 

-V.(E  xH)+H  •  V  x  E  =  E  •  f  +  E^_ 

Substituting  the  value  of  V  x  E  from  Faraday's  law  (2A)  and 
regrouping  the  terms: 

-V*(E  xH)  =  E  •/+£  •  +  ft  »  ^L 

By  integrating  over  the  volume  of  the  system  and  transforming  the 
volume  integral  on  the  left  side  into  a  closed  surface  integral  covering  the 
volume  (using  Gauss'  divergence  theorem),  the  relation  (5.1)  is  obtained. 

There  is  a  subtle  difficulty  in  using  (5.1)  to  the  case  when  the 
accelerated  body  (projectile)  is  moving.  We  assume  that  the  surface  of 
integration  is  driven  together  with  the  projectile,  also  e  and  p  remain 
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"substantially"  constant  (or  in  other  words,  a  point  in  movement  maintains 
its  properties3).  Then  the  relation  (5.1)  remains  valid  with  the  assumption 
that  the  partial  derivatives  of  D  and  B  with  respect  to  time,  are  replaced  by 
their  flux  derivatives  with  respect  to  time 


3D 

dt  dt 

(5.4) 

dB  AA 

dt  ^  dt 

(5.5) 

where,  the  flux  derivatives  of  the  fields  d  and  B  are  [25] 

d.D  pjr\  -  ,  -  s, 

~fjj~  ~  ft  +  u  div  D  +  curl  (D  x  u) 

(5.6) 

d .  B  35  _  ,  _  v 

-  ft  +  u  div  B  +  curl  (B  x  u  ) 

(5.7) 

and  using  the  other  two  Maxwell’s  equations: 

div  ^  =  A 

(5.8) 

and 

div  B  =  0 

(5.9) 

d ,  D  3ft  ,  _  v 

~d7~  ~  dt  +Am+Vx(Dxm) 

(5.10) 

3  If  p  remains  "substantially"  constant 
velocity  of  the  point. 


=  ~-  +  ugrad/i 
dt  dt 


where  u  's  the  instantaneous 


58 


=  f-  +  Vx(B 


X  u  ) 


(5.11) 


5.2  GLOBAL  (GENERALIZED)  FORCES  OF  ELECTROMAGNETIC  ORIGIN 


One  important  way  in  which  mechanical  forces  of  electromagnetic 
origin  responsible  for  driving  electromagnetic  accelerators  can  be  evaluated  is 
by  employing  the  methods  of  analytical  (Lagrangian)  dynamics  consistent 
with  the  first  principle  of  thermodynamics  and  with  Maxwell's  equations. 

For  elementary  displacements  between  two  states  of  an  insulated 
thermodynamic  system,  the  first  principle  postulates  that  the  sum  of  the 
elementary  mechanical  work  performed  by  mechanical  forces  dL  and  of  the 
elementary  increase  in  the  thermal  energy  dQ  must  equal  the  change  in  the 
internal  energy  of  the  system,  d  W  . 

~d  W  =  dL  +  dQ  (c  I?) 


In  the  presence  of  the  electromagnetic  field,  the  terms  dependent  on 
state  variables  of  the  electromagnetic  field  can  be  separated  from  the  internal 
energy  of  the  system.  Repeating  the  statement  from  the  last  section,  the 
electromagnetic  energy  of  a  system  is  that  part  of  the  internal  energy  which 
varies  with  the  variation  of  the  electromagnetic  state  variables,  remaining 
constant  when  they  remain  constant. 

As  a  consequence,  the  balance  of  energy  (5.12)  can  be  written  separately 
for  the  electromagnetic  energy  in  which  case  the  elementary  mechanical  work 
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is  produced  by  mechanical  forces  of  electromagnetic  origin  and  the  thermal 
energy  accounts  for  the  joule  (and  hysteresis)  losses  in  the  system.4 

In  rewriting  the  energy  balance,  we  are  using  the  fundamental  identity 
for  electromagnetic  energy  (Al).  Then,  for  the  infinitesimal  time  dt : 


-d  W  =  P,dt  +  Psdt  +  dL 


(5.12a) 


where 


p,  =  JjJE  *  Jdv 


represents  the  joule  losses  in  the  system  (W) 
Ps  =  -  $j)(ExH  )•  ds 


represents  the  electromagnetic  power  input,  flowing  into  the  system  through 
its  boundary  surface  Sv  (the  external  normal  to  the  surface  is  by  convention 
the  positive  one). 

The  chosen  transformation  -  in  the  sense  of  Lagrangian  dynamics  - 
must  permit  an  expedient  evaluation  of  the  mechanical  work  dL  and  from 
it,  of  mechanical  forces  of  electromagnetic  origin.  Such  a  transformation  is 
defined  by  the  constraint: 


P  dt  +  Psdt  =  0 


(5.12b) 


then 


4  Since  the  electromechanical  systems  are  notoriously  efficient,  a  lossless  model  of  the  system 
gives  sufficient  accuracy. 
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-d  W 


=  dl 

pJ+p,  =o 

V 


(5.12c) 


Let  us  express  the  constraint  (5.12b)  by  using  the  identity  (5.1): 


P  +  P 
I  s 


-d.D  _  d.D 
E  —  +  H  -L 


dt 


dt 


d  v  =  0 


(5.1 2d) 


In  order  to  impose  the  constraint  the  integrand  in  the  left  side  of  the 
equation  (5.12d)  must  be  zero.  The  two  instances  in  which  this  happens  are 
(a)  the  fields  (electric  and  magnetic)  are  zero,  instance  with  no  significance 
and  (b)  the  flux  derivatives  are  kept  zero  during  the  transformation 
(interesting!). 


dfD 

dt 


-  =  0  leads  to  the  constraint  Vti .,  ~  const 


dfB 

—‘z — =0,  ,  <Pm  =  const . 

dt  leads  to  the  constraint 


Substituting  (5.13a  and  b)  in  (5.12c) 


-  d  W 


=  dL 


<D_  =c onst  . 


(5.13a) 

(5.13b) 


(5.12e) 


actually  splitting  the  energy  dW  in  its  electric  and  magnetic  counterparts: 


dL=  -  dWe 


-  dWn 


(5. 120 


For  an  electromechanical  system  with  n  degrees  of  freedom 

h 

d  L  =  X  Fkd  xk 

k  =1 


(5.14) 


I 
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where  Fk  xk  are  the  generalized  Lagrangian  force  and  generalized 
Lagrar.gian  coordinates,  respectively. 

By  varying  only  the  x,  coordinate: 
dx;  *  0  dxk  =  0  (k  *  j) 

(dWm 

{ *, 

=ccnst 

In  summary,  the  steps  followed  in  order  to  find  the  global,  generalized 
Lagrangian  force  acting  along  the  j  th  generalized  coordinate  are: 

a.  The  system  is  constrained  to  change  such  as  [P.  +  Pc  ]  =  0; 

b.  Among  all  possible  virtual  movements,  the  variation  of  onlv  one 
generalized  coordinate  is  chosen; 

c.  The  expression  found  for  the  force  is  then  valid  independently  of 
the  variation  of  fluxes  or  currents  (T'  .e  constraints  of  constant 
electric  and  magnetic  flux  densities  describe  the  transformation 
employed  in  order  to  derive  the  expression  for  the  force  and  does 
not  affect  the  application  of  this  expression  in  the  situation  of 
variable  flux,  for  example.) 

d.  Choosing  the  right  generalized  coordinates  (their  number  being 
fixed  at  the  number  of  degrees  of  freedom  of  the  system)  can 
provide  a  special  insight  into  the  electromagnetic  accelerator 
operation. 


Some  pairs,  generalized  coordinate  -  generalized  force  are  shown  in 


table  A.I. 


TABLE  A.I. 

Generalized  Coordinate 

Generalized  Force 

Displacement,  dx  (m) 

Force,  F  (N) 

Rotation  angle  d  a  (rad) 

Torque  T  (N*nA 

Surface  ds  s  (m2) 

Superficial  Tension  t(N/m) 

1 

Volume  d  v  (m3) 

Pressure  p  (N/m2) 

i 

For  example,  by  expressing  the  energy  as  a  function  of  thp  volume,  the 
generalized  force  results  as  the  (magnetic)  pressure  on  the  accelerator  system. 

Expressing  the  magnetic  energy  as  functions  of  electric  flux  and 
magnetic  flux  (actually,  magnetic  flux  linkage)  as  the  independent  variable 
may  not  be  the  most  convenient  choice.  Many  times  the  use  of  current  as 
independent  variable  in  the  energy  expression  represents  a  more  natural  and 
convenient  alternative.  The  Legencre  transformation  from  analytical 
mechanics  is  the  tool  changing  an  energy  function  of  a  given  set  of  variables 
into  a  new  energy  function  of  a  new  set  of  variables,  the  old  and  the  new 
variables  being  related  to  each  other  by  a  joint  transformation.  The 
transformation  has  the  remarkable  property  that  is  entirely  symmetrical  on 
both  systems  (the  same  transformation  that  leads  from  the  old  to  the  new 
system  leads  back  from  tne  new  to  the  old  system)  Since,  in  this  treatment. 
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only  magnetic  acceleration  will  be  considered,  the  new  energy  function  is  the 
magnetic  coenergy  N  m  introduced  through  Legendre's  transformation: 


W  (I  tn)  =  0mI  -  W  (0m,xyx2/....xn) 


(5.16) 


Then,  the  energy  balance  becomes: 


dL  =  dW 


(5.17) 


and  the  generalized  force  along  the  xj  coordinate: 


dW'(I  ,x1,x2,...xn 
dx 


(5.18) 


5.3  LOCAL  (BODY)  FORCES  OF  ELECTROMAGNETIC  ORIGIN  AND  THE 
FLUX  DERIVATIVE  METHOD 

Electromagnetic  accelerators  and,  especially,  the  projectile  operate  at 
extreme  stress  levels.  For  the  same  global  force,  a  design  can  fail  or 
successfully  survive  depending  on  how  judicious  the  body  forces  are 
distributed. 


The  same  methodology  which  led  to  finding  the  generalized  force, 
starting  from  Maxwell's  equations  and  being  consistent  with  the  first 
principle,  can  be  also  used  to  find  the  electromagnetic  body  forces. 

A  transformation,  in  the  sense  of  Lagrangian  dynamics,  defined  by  the 
constraint  (5.12b)  was  chosen  in  the  form  of  equation  (5.12e): 


-  dW' 


(5.12e) 


104 


In  terms  of  power  (differential  energy  d  W  in  incremental  t’me  dt  )■ 


dt 


V  =amst 

t 

0  =const  . 


dJL 

dt 


(5.19) 


Where  f  is  the  body  force 


&) 


in  the  volume  d  v  element  and  u  is 


m 


the  velocity  of  the  element  — .  The  energy  tV  ( /  )  can  be  substituted  by  the 


volume  integral  of  the  energy  density  w  [  3 


d_ 

dt 


|  \ \\zvdv  =  JJ  J/  *udv 


*P  =c onst  . 

0  *=COTXSt 


(5.19a) 


In  order  to  find  the  expression  for  the  body  force  (5.19a)  must  be 
transformed  in  a  way  to  obtain  the  vectorial  factor  which  multiplies  the 
velocity  u  in  the  right  hand  member  of  the  equation.  Splitting  the  energy 
density  into  its  electric  and  magnetic  counterparts  yields: 


JJ \f*udv 


chve 


W[(* 


)  (50  ' 

ao  ^  iS_m  o 


a  v 


(5.19b) 


(We  have  a  derivative  of  a  volume  integral  on  a  domain  with  moving 
bodies.  An  assumption  is  made  that  through  the  surface  which  surrounds 
the  volume  of  the  integration  domain  bodies  do  not  "enter"  or  "come  out"  - 
or,  in  other  words,  the  surrounding  surface  is  being  "driven"  by  polarized  or 
magnetized  bodies.) 
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The  energy  densities  in  5.19b  are  given  by  5.1b  and  5.1c.  Taking  the 
derivative  by  parts 


jjjf'udv  =- J  J  J 


E1# 


it 


e2  4.  a  (SB 

Tw +  H  br. 


f£=, o 

rft 


H_^/i 
2  oT 


u  t; 


(5.19c) 


Observation:  Considering  that  the  material  properties  e  and  p  remain 
substantially  constant,  or  in  other  words,  a  point  in  movement  main,. .ns  its 
properties,  yields: 


=  ff  +  “  *  Srad£  =°  °r  grade;  ~=-u*gradp 


(5.12) 


The  condition  of  constant  flux  (electric  and  magnetic)  imposes  the 
constraint  of  zero  flux  derivatives  for  D  and  B  : 


d ,  D  ar>  _  _  ,  _  , 

^  +  m  div  D  +  curl(D  x  u  )  =  0 


(5.6) 


and 


_  dB_ 
dt  dt 


+  u  div  B  +  curl(B  +  u)  =  0 


(5.7) 


and,  by  using  Maxwell's  equations: 
dir;  D  =  ft  and  div  B  =0; 

=-[cur/(D  x  u)+  pu] 


d  D 

_l _ 

dt 


(5.6a) 
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f-)  =  -[curl(S  x  u)] 


(5.7a) 


Substituting  (5.6a)  and  (5.7a)  and  (5.19d)  back  into  (5.19c): 
JJJ/  *  “  =  [$,£  •  u  +  E  •  curl(D  x  u)  - 

V  V 

[~Srad£)»  u  +  H  •  curlitS  x  u) 


(5.20a) 


In  the  right  side,  two  terms  still  do  not  have  u  as  an  explicit  vectorial 


factor.  In  order  to  transform  them  we  use  the  vectorial  identities: 
div[E  x  (D  x  «)]  *  £  «  curl(D  x  w  )  -  (D  x  u)  •  curl  E 

div[H  x  (B  x  u)j  =  H  •  curl(D  x  u)  -  (B  x  u)  •  curl  H 


(5.20b) 


(5.20c) 


Also,  in  the  condition  of  constant  electric  and  magnetic  flux,  the 
Faraday's  law  and  the  magnetic  circuit  law  become 


curl  E  =  — —  =0;  curl  H  =  J  +  =  J 


(5.20d  and  e) 


The  substitution  of  equations  (5.20a,  b,  c,  d,  and  e)  into  (5.20-a)  leads  to  a 
remarkable  expression: 

11  \f  •  =1  jl[pE  -  —grade  +  (/  x  B )  -  ^grad\i  •  udv 

-|)|>[£  x(d  x  n)  +  H  x(B  x  n)]«ds 

’•  (5.200 
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The  surface  integral,  which  results  from  the  two  divergence  terms 
(5.20b)  and  (5.20c)  becomes  zero  after  applying  Gauss'  theorem  since  we  chose 
the  surface  surrounding  the  domain  such  that,  u  =  0  on  it.  Then,  the  body 
force  f  : 


f 


grade  +  /  x 


gradu 


(5.21) 


where  the  electric  field  intensity  E  \m  )  and  the  magnetic  flux  density 

B  (T  )  can  be  seen  as  volume  force  densities  exerted  by  the  field  on  the  unit 

C_ 

element  of  free  charge  A  m3  and  respectively  of  moving  charge,  / 


A 

(conduction  current  density  m2  •) 


The  other  two  terms  proportional  to  the  gradients  of  electric 
permittivity  and  magnetic  permeability  represent  the  body  forces  exerted  by 
the  electric  field  on  the  polarizable  materials  and  respectively  by  the  magnetic 
field  on  magnetizable  materials. 


In  summary,  the  steps  followed  in  order  to  find  the  body  force  (the 
electromagnetic  force  density  per  unit  volume)  are: 

a)  The  system  is  constrained  to  change  such  as  the  electric  and 
magnetic  flux  remain  constant  during  the  transformation 
consistent  with  the  first  principle  of  thermodynamics  and 
Maxwell's  equations; 
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b)  The  constancy  of  the  electric  and  magnetic  flux  is  imposed  by 
equating  to  zero  the  flux  derivatives  of  the  field  vectors  D  and 
B  ; 

c)  The  body  force  is  obtained  as  a  vectorial  factor,  multiplying  the 
local  velocity  u  in  the  expression  for  electromagnetic  power; 

d)  The  body  forces,  thus  obtained  are  universally  valid  for 
macroscopic,  nonrelativistic  electrodynamics; 

3)  In  the  above  derivation,  the  variation  of  e  and  M  with  respect  to  the 
volume  density  of  the  materials  (T)has  been  neglected.  If  such 
variation  exists,  the  electro-  and,  respectively,  the  magneto-striction 
terms  will  appear  in  (5.21). 


5.4  THE  ELECTROMAGNETIC  STRESS  TENSOR  AND 
ELECTROMAGNETIC  IMPULSE 


The  Maxwell's  electromagnetic  stress  tensor  was  actually  based  on 
Faraday's  idea  of  an  intermediary  medium  with  special  electric  properties  for 
transmitting  the  force.  Let  us  again  assume  the  domain  /  ^  surrounded  by 
the  surface  Sv  (Figure  4)  The  force  Fz  =  -Fn  is  introduced  to  maintain  the 


equilibrium  when  the  imaginary  elastic  strings  are  cut.  This  way,  the  global 
electromagnetic  force  F  (N  )  is  equivalent  to  the  stresses  ^"(^2)  integrated 
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over  the  area  ,  which  in  turn,  balances  the  body  forces  of  electromagnetic 
origin  /  (~i)  integrated  over  the  volume  / 

F  =  =$>$TJs 

v  (5.23) 

The  replacement  of  body  forces  by  stresses  (using  the  Gauss’ 
transformation  of  a  volume  integral  into  a  surface  integral)  is  performed  for 

SD 

the  particular  case  of  a  stationary  regime  -  0  and  absence  of  material 

motion  «  =  0  -5  At  the  end  of  this  paragraph,  the  nonstationary  regime  will 
be  considered  while  evaluating  the  electromagnetic  impulse  (so  necessary  in 
calculating  different  aspects  of  the  recoil  in  electromagnetic  guns). 

Starting  from  equation  (5.21)  we  will  try  to  express  the  body  force  of 
electromagnetic  origin  as  a  divergence,  in  order  to  find  an  equivalent  stress 
distribution 

_  p2  _  _  pi2 

f  =  fl,E  -  —grade  +  J  x  B  -  —  grad^i  = 

EdivD  -  grad  +  £grad  +  curlH  x  B 

-grad 

in  which  the  last  two  terms,  added  for  symmetry,  have  zero  value. 


H  ______ 

+  figrad  —  +  curJE  x  D  +  HdivB 


(5.21b) 


5  The  expression  for  the  body  force  (A-21 )  was  obtained  in  the  previous  paragraph  under  the 
constraint  of  zero  flux  derivatives  which  in  case  of  absence  of  motion  w  =  0  yields  the  equations 
curlH  =  J  and  curlE  =  0 . 
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By  using  the  vectorial  identity: 
(V  x  E  )  x  D  = 

-eE  x  (V  x  E  )  =  -e[v(-y)  ~  V)E  ] 


and 

(V  x  H)  x  B  = 

-MH  x  (V  x  B  )  =  -  (H  V)H 


(5.24a) 


(5.24b) 


the  volume  force  density  expression  becomes: 

/  =  /«.+ w= 

£rfc;D  -grad  (y-)  +  (y-)-  egrad  (y)  +  (Dgrad  )E 

+  Hd  ivB  -  grad  (-y^~]  +  (y-)  -  /igrarf  (y~)  +  ( Bgrad  )H  ^ 

In  relation  (5.21c)  there  is  a  perfect  symmetry  between  the  electric  and 
magnetic  field  quantities.  The  right  side  of  the  equation  (5.21c)  expresses  the 
tensoral  divergence  of  a  dvdic: 

The  x  components  of  the  electric  and  magnetic  body  forces  in  this 
relation  are: 

ftlx  =  Exd  ivD  +  D  •  gradEx  -  ^-(y-)  =  div(Ex  x  D  )  -  div^f  y-j 


(5.21d) 


Ill 


fmag*  =  HzdivB  +  B  •  gradHx  -  (^y“)  =  div(Hx  x  B  )  -  div  i  ■— 


(  _  ^rr2 
~2 
(5.21e) 


considering  also  the  y  and  2  components,  we  can  express,  finally,  the  body 
force  as  a  tensorial  divergence: 


/  =  L.+  Uag  =  1  div  Tx+ ]  div  Ty  +  k  div  Tz 


(5.25) 


where  f ,/  ,  k  are  the  unit  vectors  of  the  three  axes  of  coordinates  x  ,  y  ,  2 


and 


rT  = 


r  = 

x  y 


T,  = 


e,d  +  h,b  -  f  (^-  +  ^ 


EyD  +  HyB  -  f 


2  2 


'g£a  MH2 
V  2  2 


Ez  D  +  B  -  fc  ( ^ 


(5.25a) 


(5.25b) 


(5.25c) 


integrating  by  components 

F=JJJ^n=JJJf  din  Txd  v  +  J  J  J  f  div  Tydv  +  j  \ \k  div  Tzdv 

V  V  v  V 

(5.25d) 


and  using  Gauss'  theorem  we  can  transform  the  volume  integral  of  a 
divergence  into  a  flux  through  a  closed  surface  integral: 

$$[*"  (jx  •  n  )  +  f  (fy  •  n  )  +  k  (fz  •  n  )]ds  =  |>jxf„ds 


(5.26) 


112 


where 


T  = 


(T  T  T  ^ 

I  1  it  1  xy1  xz 

TyxTyyTyz 

VT«  V. 


is  the  electromagnetic  stress  tensor 
and 

fn=T  *  n  =  (i~;fz  +  f;  fy  +  k;Tz)  ♦  n 
F  =  <j)|)r  n  f  dv  =  j  j  J  f  dv 

6  V  V 

v 


Then,  the  normal  component  of  the  stress  tensor  is: 
f„=  E  (d  •  w  )  +  H  (b  •  n)  -  n(zvel  +  wm) 


and,  the  expression  for  the  Maxwell's  tensor: 


T  = 


(EZDX  +  HZBX  -  zu  ExDy  +  HXB} 
EVDX  +  HyBx 


EXDZ  +  HXB, 


y  x  •  **  y  ^x 

Ez  Dx  +  Hz  Bx 


EyDy  +  HyBy  ~  W  EyDj  +  HyBj 

Ex  D  y  +  Hz  B  v 


EzDz  +  HZBZ  -  zv 


where  w  -  is  the  electromagnetic  energy  density 

are  interested  only  in  the  magnetic  terms 


ify  (i5)- 


fm=J  x  B  -  ~ gradfj.  =  div  T„ 


(5.23a) 


(5.27) 


Since  we 


(5.27a) 


then 
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Tm  •  n  =  H  (B  •  *  B  ) 


(5.27b) 


where  n  is  the  positive  normal  to  the  surface.  As  a  rule,  H  the  intensity  of 
the  magnetic  field  bisects  the  angle  between  n  and  f 

Electromagnetic  Impulse 

Suppose  that  the  regime  is  not  stationary  as  was  assumed  at  the 

_  _  3rS 

beginning  of  the  paragraph  (5.20e).  Then  curlH  =  J  +  .  In  this  case,  the 

equation  (5.23)  must  be  modified: 

j  j  Sfdv  -  =-jy-JJJ(D  x  B)dv 

v  «.  0  (5.29) 


or 

F  +  J  J  J  (D  x  B  )d  v  =  |>|>( -T  )(-n)ds 
0  *. 


(5.29a) 


If  we  consider  an  isolated  system,  which  extends  to  infinity,  then  the 
surface  integral  disappears  since  ^  goes  to  zero  faster  than  ds  is  increasing 
for  an  R  —>°o.  This  yields: 

F  +jfl!S(v  x  B)dv  =0 


(5.29b) 


114 


The  resultant  force  in  an  isolated  mechanical  system  interacting  with 
an  electromagnetic  system  is  zero;  expressed  in  terms  of  the  mechanical  and 
electromagnetic6  impulse: 


mech  . 


+ 


1JJ(S 


(5.30) 


or  the  sum  of  the  impulse  is  constant 


+  x  B)dv 

v 


=  const. 


In  general. 


-t  )(-« w* 8  <5-30a) 

, 

The  expression  (5.30a)  is  important  in  evaluating  the  design  questions 
related  to  recoil  in  electromagnetic  guns.  The  volume  integral  of  the 
electromagnetic  force  density,  the  Maxwell's  stress  tensor,  and  the 
electromagnetic  impulse  become  tools  in  defining  and  calculating  globally 
and  locally  the  recoil  in  electromagnetic  accelerators. 


6  The  density  of  electromagnetic  impulse  is  D  x  B  =  £flE  x  H  =  £[JS 
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CHAPTER  6 

SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 
6.1  SUMMARY 

It  is  considered,  in  this  work,  that  one  of  the  reasons  why  the 
electromagnetic  hypervelocity  launchers  and  associated  power  supply 
technologies  have  not  reached  their  expectations  is  the  fact  that  most  of  the 
activity  in  this  domain  has  been  focused  on  simple  railgun  accelerators  in 
conjunction  with  the  questionable  choice  for  homopolar  generators  known 
for  their  inherent  low  voltage  and  low  power  density.  In  order  to  increase  the 
power  densities,  the  homopolar  machines  have  to  be  used  in  conjunction 
with  heavy  and  voluminous  inductors  in  which  the  magnetic  energy  is 
transiently  stored  and  explosive  opening  switches  which  actually  provide  the 
high  voltages  necessary  for  driving  high  currents  into  railguns. 

In  order  to  drastically  improve  the  performance  of  hypervelocity 
electromagnetic  launchers  the  solution  is  to  use  as  power  supplies  complex 
rotating  flux  compressors  (beyond  so-called  compulsators)  in  conjunction 
with  heteropolar  electromagnetic  launchers  with  variable  structures. 

The  computation  effort  for  the  design  of  such  complex  magnetic  flux 
compressors  and  heteropolar  hyper  velocity  accelerators  is  tremendous  in 
order  to  determine  the  time  evolution  of  complex  interacting  electromagnetic 
fields. 
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For  the  demanding  analysis  and  design  of  flux  compression  structures, 
a  three-dimensional,  finite  element  method  formulation  for  the 
electromagnetic  field  diffusion  is  presented. 

The  formulation  is  based  on  the  Galerkin  technique  since  in  the 
domain  of  parabolic  partial  derivative  equations  no  variation  principle  is 
possible.  This  formulation  requires  large  computational  efforts  since  it  solves 
the  transient  problem  for  all  three  components  of  the  vector  potential  (as 
compared  with  axisymmetric  f.e.m.  codes  solving  for  one  component  only). 

One  of  the  original  features  of  the  formulation  is  the  introduction  of 
superficial  currents  (sheets  of  currents)  with  linear  current  density,  K , 
flowing  through  the  faces  of  each  volume  element  (brick,  prism,  etc.)  and 
compensating  for  the  discrete  nature  of  the  finite  element  method  and  for 
other  nonuniformities  and  discrepancies  such  as  in  prescribing  the  forcing 
functions  (current  densities,  for  example)  element  by  element  enforcing 
point-wise  the  solenoidality  of  flux  density,  B,  current  density,/,  and 
magnetic  vector  potential,  MdivB  =0;  divj  =  0;  divA  =0). 

The  superficial  currents  with  linear  current  densities  K  are  introduced 
in  the  f.e.m.  formulation  as  Lagrange  multipliers,  solenoidality  of  the  three 
vectors,  8,  /,  and  A ,  being  the  constraint  imposed  by  the  Lagrange 
multipliers. 

The  unfolding  in  time  of  the  flux  compression  phenomena  occurring 
in  a  three-dimensional  flux  compression  structure  (pulsed  coil  in  front  of 
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consecutively,  producing  remarkable  velocities.  Armature  reaction  will 
introduce,  however,  complex  three-dimensional  field  distributions. 

In  the  domain  of  the  power  supplies,  complex  rotary  flux  compressors 
as  those  proposed  by  Driga  and  Fair  in  [35]  must  be  analyzed  by  the  analytical 
and  numerical  methods  treated  in  the  present  work. 

One  of  the  features  desired  to  be  incorporated  into  the  f.e.m.  3-D  code  is 
the  relative  motion  between  the  rotor  and  stator  for  the  flux  compressors  and 
between  the  barrel  and  the  projectile  for  the  accelerators.  In  order  to 
accommodate  such  problems,  the  f.e.m.  computer  code  must  be  adaptive  with 
possibility  of  continuous  change  in  e  structure  of  formulation. 

The  program  must  have  an  "initial"  mesh  changes  according  to  the 
movement  of  the  projectile. 

The  algorithm  must  assume  two  rigid  bodies,  both  generating  magnetic 
fields  moving  with  relation  to  one  another  in  an  accelerated  or  decelerated 
manner. 

A  mesh  of  brick  elements  is  generated  around  each  bo  ’y  and  those 
meshes  meet  along  a  mesh  interface,  I.  To  provide  for  a  smooth  sliding  of 
one  element  with  respect  to  another,  a  special  five  mode  interface  element  is 
used  [36].  The  associated  shape  functions,  ^(C/7?),  are  defined  on  a  master 
element,  £}.  The  moving  local  coordinates  on  a  typical  element  are  then 
given  by  an  isoparametric  map. 


118 


While  codes  with  adaptive  mesh  exist  for  solving  problems  in  other 
domains  of  technology  (such  as  supersonic  flow  in  domains  with  moving 
bodies),  they  do  not  exist  for  electromechanical  systems  where  they  will  have 
an  important  impact  on  analyzing  highly  transient  systems. 
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